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ABSTRACT 

A large sample of spectroscopically confirmed star-forming galaxies at redshifts 1.4 < z spcc < 3.7, 
with complementary imaging in the near- and mid-IR from the ground and from the Hubble and Spitzer 
Space Telescopes, is used to infer the average star formation histories (SFHs) of typical galaxies from 
z ~ 7 to 2. For a subset of 302 galaxies at 1.5 < z spoc < 2.6, we perform a detailed comparison of 
star formation rates (SFRs) determined from SED modeling (SFRs[SED]) and those calculated from 
deep Keck UV and Spitzer /MIPS 24 /j,m imaging (SFRs[IR+UV]). Exponentially declining SFHs yield 
SFRs[SED] that are 5 — lOx lower on average than SFRs[IR+UV], indicating that declining SFHs may 
not be accurate for typical galaxies at z > 2. The SFRs of z ~ 2 — 3 galaxies are directly proportional 
to their stellar masses (M*), with unity slope — a result that is confirmed with Spitzer/IRAC stacks 
of 1179 UV-faint (K > 25.5) galaxies— for M * > 5 x 1O 8 M and SFRs > 2M yr- 1 . We interpret 
this result in the context of several systematic biases that can affect determinations of the SFR-Af* 
relation. The average specific SFRs at z ~ 2 — 3 are remarkably similar within a factor of two to 
those measured at z > 4, implying that the average SFH is one where SFRs increase with time. A 
consequence of these rising SFHs is that (a) a substantial fraction of UV-bright z ~ 2 — 3 galaxies 
had faint sub-L* progenitors at z > 4; and (b) gas masses must increase with time from z — 7 to 2, 
over which time the net cold gas accretion rate — as inferred from the specific SFR and the Kennicutt- 
Schmidt relation — is ~ 2 — 3x larger than the SFR . However, if we evolve to higher redshift the SFHs 
and masses of the halos that are expected to host L* galaxies at z ~ 2, then we find that < 10% 
of the baryons accreted onto typical halos at z > 4 actually contribute to star formation at those 
epochs. These results highlight the relative inefficiency of star formation even at early cosmic times 
when galaxies were first assembling. 

Subject headings: dust, extinction — galaxies: evolution — galaxies: formation — galaxies: high- 
redshift — galaxies: star formation 



1. INTRODUCTION 

In the last few years, it has become a standard prac- 
tice to decipher the physical characteristics of distant 
galaxies by fitting broadband photometry with spectral 
synthesis models. Stellar population modeling, as it is 
called, has been aided by the availability of deep imag- 
ing in extragalactic fields across a large baseline in wave- 
length. Comparison of the broadband spectral energy 
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distribution (SED) of galaxies with that of a population 
of stars with a given initial mass function, star formation 
history, age, dust reddening, and metallicity, can there- 
fore yield important insights into the physical proper- 
ties of high-redshift galaxies. This modeling has become 
more sophisticated, with some versions allowing for the 
presence of strong emission lines (or simultaneously fit- 
ting for such l ines) that may affect the broadba nd pho- 
tometry (e.g.. ISchaerer fc de Barrosl I2009L [2010h . Other 
models incorporate the full stellar and dust SEDs in or- 
der to derive self-consistently the reddening of starlight 
based upon direct dust indicators (e.g., such as the mid- 
or far- infrared dust continuum), thus accounting for the 
total energy budget when fitting for the stel l ar popula- 
tions (jGordon et alJl200H iMlsselt et alj|200lt iNoll et all 
2009). The latter have somewhat limited use for high- 
redshift galaxies since it is only for the most infrared 
luminous and dusty galaxies at z > 2 that individual 
detections at mid and far-infrared wavelengths are at- 
tainable, thus allowing the modeling of the full IR SED. 

While there has been much progress in developing ever- 
sophisticated methods of fitting the stellar populations 
of distant galaxies, the one fundamental obstacle that 
affects most of these methods is the inherent degeneracy 
between the star formation history, age, and dust red- 
dening, even when the redshift of the galaxy is known 
beforehand (e.g., from spectroscopy). Lack of redshift 
information will of course only further hinder one's abil- 
ity to robustly determine these quantities. It is difficult, 
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if not impossible, to reliably disentangle these effects 
based on broadband photometry alone, even with the 
deepest optical and near-IR data, as has been discussed 
in the first investigations that mo deled the stellar pop- 
ulations of high-redshift galaxies (iSawicki fc Yed 119981 : 
iPapovich etaLll200H iShaplev et al.ll2001[ ). Full spectral 
energy distribution modeling of the stellar and dust com- 
ponents can break some of this degeneracy, but can also 
add a new layer of complication given the increase in 
number of free parameters that describe the stellar pop- 
ulation and the dust properties and the spatial distribu- 
tion of that dust with respect to the stars in a galaxy. 

Finally, there are some inherent uncertainties in SED 
modeling that will likely never be fully resolved. In par- 
ticular, even in the best case with deep UV through 
near-IR photometry, the data are still insufficient to dis- 
tinguish simple star formation histories (such as those 
parameterized as monotonic exponentially declining, ris- 
ing, or constant functions) from more complicated ones 
that include multiple generations of bursts. In contrast 
with fossil studies of near by resolved stellar populations 
(e.g. JWilliams et al.l2010f) . and in the absence of detailed 
spectroscopic abundance measurements, it is difficult to 
work "backwards" from the integrated light of the stellar 
populations in a galaxy to a unique set of star forma- 
tion histories for that galaxy. Nonetheless, simple star 
formation histories that vary monotonically with time 
have been widely used to infer the stellar population pa- 
rameters for high-redshift galaxies. The most commonly 
adopted function is one in which the star formation rate 
(SFR) of a galaxy declines exponentially with time, as 
would be predicted from a closed box model of galaxy 
evolution dSchmidtl H95l iTinslevI H980ft . Such exponen- 
tially declining models have been popular as they repro- 
duce the optical /nea r-IR colors of local spiral galaxies 
(jBell fc de JoneKOOOl ) and appear to reproduce the over- 
all evolution in the SF R density at redshifts z < 2 (e.g., 
iNaeamine et al.l feOOOft . As surveys of dropout selected 
samples push to increasing redshifts, however, it has be- 
come clear that galaxies at z > 2 have SFRs and stellar 
masses that are inconsistent with them having formed 
stars according to an exponentially declining or constant 
star formation (CSF) history prior to the epoch during 
which they are observed. 

Recently, there has been substantial interest in the pos- 
sibility that high-redshift galaxies in general may fol- 
low "rising" star formation histories, where the SFR 
increases exponentially or linearly with time. Circum- 
stantial evidence for such rising star formation histories 
comes from pred ictions of cosmological hydrodynamic 
simul ations (e.g., iFinlator et al.l 1201 It IWeinman n et al.l 
120111 1 . the presence of SFR v ersus stellar mass correla- 
tions at redshifts z > 1 (e . g.. [Renzinil [20091 : iStark et al.1 
120091: iGonzalez et al1l2010t ILee et al.ll2011a| ) , and the in- 
crease in th e SFR density per cq moving volume at early 
times (e.g.. IPapovich et aLll2011|l . 

While such simple monotonic functions are unlikely to 
capture the full diversity and complexity in the star for- 
mation histories of galaxies, we can still make progress 
by addressing the average statistical properties of galax- 
ies across a wide range in redshift, or lookback times, 
to effectively look back into the history of star forma- 
tion and thus attempt to deduce the way in which aver- 
age galaxies are evolving. One method is to use cluster- 



ing measurements and halo abundance matching to infer 
a "d uty cycle" for star formation on a stati stical basis 
(e.g.. lAdelberger et alHi998t ILee etldl 12009ft . Another 
method is to use multi-wavelength indicators of redden- 
ing and total SFR to constrain certain parts of parameter 
space spanned by the SED fitting parameters. This ap- 
proach is simpler than that taken by studies that treat 
the dust properties and distribution of dust with respect 
to the stars as additional free parameters in modeling 
the full spectral and dust SEDs. Further, the advantage 
of performing direct comparisons between independent 
indicators of SFR and those derived from SED fitting is 
that the method can be applied to individual galaxies, 
so long as they are detected (or have meaningful upper 
limits) at mid-infrared wavelengths, in order to indepen- 
dently measure the fraction of dust obscured light. And, 
unlike stellar masses which are typically exclusively mea- 
sured from the rest-frame near-IR light (which can also 
have a significant contribution from current star forma- 
tion), there are many independent methods of estimat- 
ing SFRs from continuum emission (UV, infrared, radio) 
or nebular line emission (e.g., Ha, Paa), thus allowing 
one to investigate the systematics and cross check results 
from different methods. 

In this paper we investigate the typical star forma- 
tion histories of spectroscopically confirmed UV selected 
star-forming galaxies at redshifts 1.4 < z < 3.7. We in- 
corporate in our analysis deep Spitzer/MlPS 24 /im data 
that exist for a subset of galaxies in our sample, in the 
redshift range 1.5 < z < 2.6; these mid-IR data are used 
to place independent constraints on the SFRs and dust 
reddening of galaxies in our sample, quantities that are 
then compared to those obtained from the SED fitting 
given various assumptions of the star formation history. 
We then proceed to discuss this comparison in light of 
recent results at higher redshifts (z > 3) to form a con- 
sistent picture for the typical star formation history of 
galaxies during the first ~ 3 billion years of cosmic time. 
Our sample and analysis lend themselves uniquely to ad- 
dressing these broad questions because of the large num- 
ber (N — 1951) of spectroscopic redshifts in the range 
1-4 < z spcc < 3.7; the deep UV, optical, near-IR, and 
IRAC data necessary to model the stellar populations; 
the deep MIPS 24 (im data, used as an independent probe 
of dust attenuation and bolometric SFR; and the careful 
consideration of the typical assumptions in SED model- 
ing, and biases in determining the relationship between 
SFR and stellar mass. 

The outline of this paper is as follows. In Sec- 
tion [2] we briefly describe the color criteria used to se- 
lect the sample of 1.4 < z < 3.7 galaxies, and sum- 
marize the fields targeted. In addition, we present de- 
tails of the multi-wavelength data in our fields, includ- 
ing ground-based near-IR and Hubble/WFC3 imaging, 
Spitzer/IRAC imaging, and Spitzer/MlPS 24 /im obser- 
vations. The rest-frame UV through near-IR photometry 
is used to constrain the stellar populations of spectro- 
scopically confirmed galaxies in our sample, as discussed 
in Section [3] A detailed comparison between the bolo- 
metric measures of star formation obtained by combining 
the UV and MIPS 24 data, with those obtained from 
the SED modeling, is presented in Section SJ Section [5] 
focuses on modeling the "younger" galaxies in our sam- 
ple by taking into account dynamical time constraints on 
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the ages and a systematic steepening of the UV attenua- 
tion curve with younger stellar population age. The sys- 
tematic variations and random uncertainties in the ages 
and masses of galaxies in the spectroscopically confirmed 
sample are discussed in Section|6] In addition, we present 
our determination of the SFR versus stellar mass correla- 
tion, and show how Malmquist bias can affect inferences 
of the slope of this relation at high-redshift. In Section[3 
we examine the mass-to-light (M/L) ratios of z ~ 2 — 3 
galaxies at rest-frame UV through near-IR wavelengths, 
and present a stacking analysis of the IRAC data on UV 
faint galaxies lying below our spectroscopic limit. The 
M/L ratios and stacking results are then used to infer 
the stellar masses of UV faint galaxies. In Section [8] we 
discuss the implications of our results for the typical star 
formation histories of high-redshift galaxies; the progen- 
itors of z ~ 2 — 3 galaxies; and the time evolution of 
cold gas mass and net gas accretion rate with redshift. 
For ease of com parison with the literature, we assume a 
iSalpeterl (jT955) initial mass function (IMF) and adopt a 
cosmology with H — 70kms _1 Mpc -1 , il\ — 0.7, and 

n m = 0.3. 

2. SAMPLE 

2.1. Galaxy Selection and Optical Spectroscopy 

Galaxies at redshifts 1.4 < z < 3.7 were selected us- 
ing the BM, BX, and Lyman break galaxy (LBG) rest- 
UV c olor criteria (ISteidel et afl 120031 : lAdelberger et al.l 
12004 ISteidel et al.l 120041 ). The imaging data were 
obtained mostly from using the Palomar Large For- 
mat Camera (LFC) or th e Keck Low Resolution Imag- 
ing S pectrograph (LRIS; lOke et al.l 119951: ISteidel et at] 
2004). The photometry and spectroscopi c followup fo r 
this survey are d e scribe d in ISteidel et afl (|2003l 120041 ): 
lAdelberger et al.l (|2004f ). Rest-UV spectroscopy with 
Keck/LRIS was obtained for about 25% of the sample 
with 1Z < 25.5. Over all of the fields of our z ~ 2 — 3 
survey, the total numbers of photometrically selected BX 
and LBG candidates that are detected in G and 1Z with 
> 5ct significance are 25,359 and 16,655, respectively, to 
a typical depth of 1Z ~ 26.5. While most of the sub- 
sequent analysis is based on the spectroscopic sample, 
we also use the faint 1Z > 25.5 galaxy data to infer the 
stellar masses of UV faint galaxies (Section [7]). 

2.2. Near-IR Data 

Constraining the stellar population of a galaxy relies 
critically on data that bracket the rest-frame spectral 
region between ~ 3600 and 4000 A. It is at these wave- 
lengths that metal absorption lines from F, G, and K type 
stars dominate the spectrum resulting in a break around 
4000 A; an additional absorption feature at 3646 A marks 
the edge of the Balmer series and is strongest in more 
massive A stars. Both features are sensitive to age 
(though the 3646 A break reaches a maximum at inter- 
mediate ages of ~ 0.3 — 1 Gyr). To probe the strength of 
these features in z ~ 2 — 3 galaxies, we obtained J and/or 
K s imaging in 14 fields of the LBG survey, using the 
Palomar/WIRC and Magellan/PANIC instruments, to 
typical 2" aperture 3 a depths of 24.4 (K s ) and 25.0 mag 
(J). The data were reduced using IDL scripts customized 
for WIRC data, a nd photometry was per formed using 
Source Extractor (|Bertin fc Arnouta 119961) . The near- 



IR data and the red uction procedures are discussed in 
IShaplev et afl (120051) . 

In addition, we have obtained Spitzer/IRAC data for 
14 fields of the survey through General Observer (GO) 
programs in Cycles 1, 3, and 7, and through DDT time. 
When we include the GTO, IOC, and Legacy programs, 
there is a total of 27 fields in the LBG survey with IRAC 
data (Table [1]). The IRAC coverage of our galaxies typi- 
cally included either channel 1 (3.6 /mi) and 3 (5.8 /mi) or 
channel 2 (4.5 /mi) and 4 (8.0 /im), with a small fraction 
of galaxies having coverage in all four channels (e.g., such 
as galaxies in the GOODS-N field, or those that are at 
the edges of the optical images). The data were reduced 
using custom IDL scripts to correct for artifacts and flat 
field the data. Indi vidual images were mosaiced using the 
MOPEX software ijMakovoz fc Mar lead [2005h . To take 
advantage of sub-pixel dithering between individual ex- 
posures, we drizzled the final mosaics onto a grid with a 
pixel scale of 0.6" (half the native IRAC pixel scale), en- 
abling higher resolution images and finer sampling of the 
point spread function. Photometry was performed us- 
ing point-spread function (PSF) fitting to prior positions 
determined from the higher resolution optical and near- 
IR data. Errors and bias in the photometry were cal- 
culated by adding artificial sources to the IRAC images 
and recovering them using the same PSF fitting software 
used to measure photometry. The details o f the PSF fit- 
ting an d IRAC photometry are provided in lReddv et al.1 
(|2006aD . 

Finally, we have obtained 8100 sec HST/WFC3- 
F160W (//-band) imaging in 10 fields (with 14 pointings 
total) of the z ~ 2 — 3 survey, as part of the Cycle 17 GO- 
11694 program (PI: Law) . Details of the da ta acquisition 
and reduction are given in lLaw et al.l (|2012f ). Briefly, nine 
900 second exposures in each p ointing were reduced an d 
combined using MultiDrizzle (|Koekemoer et afl [2002) . 
The individual exposures are sampled onto a grid with 
a pixel scale of 0'.'08 to take advantage of the subpixel 
dithering between exposures. The typical 5 a depth of 
the combined images is ~ 27.9 AB, assuming a 0'.'2 ra- 
dius aperture. The Hubble data are particularly advanta- 
geous because the combined depth of the F160W images 
allows us to constrain the SEDs for fainter objects that 
may otherwise be undetected in the ground-based J/K s 
or Spitzer/IRAC imaging. The ground-based near-IR, 
Hubble/WFC3 F160W, and Spitzer/IRAC data, in con- 
junction with our U n G7Z optical data, are used to con- 
strain the stellar populations and stellar masses of galax- 
ies in our sample, as described in the next section. 

2.3. Mid-IR Data 

A key aspect of our analysis incorporates indepen- 
dent measurements of the SFRs of high-redshift galax- 
ies, based on direct tracers of dust emission. Six of the 
fields in our z ~ 2 — 3 survey contain deep (« 4hr) 
Spitzer/MlPS imaging at 24 /im with a typical 3 a depth 
of « 12/iJy: GOODS-N (PI: Dickinson) and West- 
phal (PI: Fazio) fields; and the Q1549, Q1623, Q1700, 
and Q2343 fields from Cycle 1 and 3 Spitzer GO pro- 
grams H These observations probe the dust sensitive 

12 There are several additional fields in our survey that contain 
Guaranteed Time Observer (GTO) MIPS imaging. These data are 
generally much shallower (typically just a few hundred seconds), 
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LBG Survey Fields with Near-IR or IRAC Data 
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a Right ascension in hours, minutes, and seconds. 
b Declination in degrees, arcminutes, and arcseconds. 
c PI: Law. 
d Pis: Steidel, Erb. 

features around rest-frame 8 /im. The MIPS data and 
reduc tion are discussed in detail in lReddv et al.l ([2006b, 
2010). Briefly, the data are flatficlded using a custom 
IDL program, and combined with the MOPEX software. 
Photometry is performed using PSF fitting to prior posi- 
tions defined by detection in the higher resolution IRAC 
data. Photometric bias and errors are estimated from 
simulations where we have added artificial sources to the 
images and recovered them using the same PSF fitting 
method. 

2.4. Subsamples and Redshift Ranges 

Throughout this paper, we use different subsamples 
of the data in different redshift ranges, with the fol- 
lowing motivations. In general, the "z ~ 2" sample 
refers to those galaxies with 1.5 < z spcc < 2.6, or 
1-4 < z S pcc < 2.7. These two different ranges arc adopted 
depending on which sample (i.e., the MIPS or SED sam- 
ple) is being used. The u z ~ 3" sample refers to those 
galaxies with 2.7 < z spcc < 3.7. Our total sample with 
available SED fits (i.e., have U n G1Z data plus at least one 
photometric point redward of the Balmer break) consists 
of 1959 galaxies with 1.4 < z spe c < 3.7. Of these, there 
are 302 galaxies with deep MIPS observations, 121 of 
which are detected individually at 24 fim with > 3 a sig- 
nificance, that allow for measurements of the rest-frame 
8 fim emission, specifically for those galaxies with spec- 
troscopic redshifts in the range 1.5 < z spcc < 2.6. This 

and are not used in this analysis. 



"MIPS" sample is used to investigate the comparison be- 
tween SED and multi- wavelength SFRs (Section [4j, and 
to investigate the relationship between SFR and specific 
SFR and stellar mass for z ~ 2 galaxies (Section I6.4|) . 
The comparison of the ages and stellar masses derived 
assuming constant and rising star formation histories is 
presented for the entire sample of 1959 galaxies in Sec- 
tion [6TTJ In Section [71 we consider the mass-to-light ra- 
tios of galaxies in our sample at different wavelengths. 
To quantify the M/L ratio at F160W and if s -band, we 
use only those galaxies with spectroscopic redshifts such 
that F160W and K s lie longward of the 4000 A break (98 
and 491 galaxies, respectively, over the redshifts ranges 
1.4 < z < 2.5 and 1.9 < z < 3.7). We use similar subsets 
of the data that have IRAC channel 1, 2, 3, or 4 data to 
quantify the M/L ratio at rest-frame 1.1 — 2.4 /im (643, 
673, 180, and 190 galaxies, respectively). The mass-to- 
light ratios at UV wavelengths are determined using 630 
and 344 galaxies with IRAC channel 1 data over the red- 
shift ranges 1.4 < z < 2.7 and 2.7 < z < 3.7, respec- 
tively. Finally, we also consider a UV-faint subsample 
with IZ > 25.5, consis ting of 1179 candidates, as dis- 
cussed in Section 17.2.21 The subsamples, their redshift 
ranges, and the number of objects, are summarized in 
Table [5J The redshift distributions of the various sam- 
ples are shown in Figure [TJ 
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Table 2 



Subsamples 


and Redshift Ranges 




Subsample Name 


Redshift Range 


ODJ 


MIPS Sample a 


1.5 < 2 ap oc < 2.6 


302 


SED Sample at z ~ 2 


1.4 < z spC c < 2.7 


1389 


sample at z ~ 


2.7 S z spC c < o.7 


570 


M/L Ratio at F160W b 


1.4 < z apC c < 2.5 


98 


M/L Ratio at A' s -band b 


1.9 < ^spcc 3.7 


491 


M/L Ratio at 3.6/im b 


1.9 < z S p C c < 3.7 


643 


M/L Ratio at 4.5 //m b 


1.9 < 2 spcc < 3.7 


673 


M/L Ratio at 5.8/jm b 


1.9 < z spcc < 3.7 


180 


M/L Ratio at 8.0(im b 


1.9 < 2 sp oc < 3.7 


190 


M/L Ratio at 1700 A with 3.6 /mi 


Coverage 1.4 < z ap ec < 3.7 


974 


Faint Sample with ft > 25.5 


BX/LBG Color Selection 


1179 



a The MIPS sample includes all galaxies that are covered by Spitzer/MIPS 24 /an 
imaging, irrespective of whether they were detected at 24 fim. 

b Includes only those galaxies in the SED sample that are detected at F160W, K s , or 
IRAC channels, and where the band lies completely redward of the 4000 A break. 
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Figure 1. Spectroscopic redshift distributions of 1959 galaxies at 
redshifts 1.4 < z sp ec < 3.7, color coded according to the subsam- 
ples considered (see Table [2}. 

In this section, we discuss the general procedure used 
to model the stellar populations of galaxies in our sam- 
ple. There are a number of assumptions that enter into 
such modeling, such as the adopted star formation his- 
tory (e.g., constant, exponentially declining, or rising), 
the imposition of a lower limit to the age of a galaxy, 
and the choice of attenuation curve. In the subsequent 
sections, we discuss and motivate our assumptions by 
utilizing the multi-wavelength data that exist in a subset 
of the fields of our survey. 

Stellar masses are inferred by modeling the broadband 
photometry of galaxies, using the full rest-frame UV 
through near-IR photometry to fit for their stellar popu- 
lations. For the fitting, we considered only those galax- 
ies that are directly detected at wavelengths longward of 
rest-frame 4000 A which, for the majority of the sam- 
ple considered here, includes the F160W, K s , and IRAC 
bands. Further, we excluded from the fitting any AGN 
that were identified with strong UV emission lines (e.g., 
Lya, CIV) or had a power law SED through the IRAC 
bands. Previous efforts to model the stell ar populations 
of gal ax ies in our sample a re described in iShaplev et al.l 
(2005), lErb et al.l (|2006bj) . IReddv et al.1 (|2006al ). and 



iReddv et al.1 ((2010) . The latest solar metallicity mod- 
els of S. Chariot & G. Br uzual (in preparati o n, her eafter 
CB11) that include the iMarigo fc Girardil (|2007t ) pre- 
scription for the thermally-pulsating Asymptotic Giant 
Branch (TP-AGB) evolution of low- and intermediate- 
mass stars are used in the fitting. The stellar masses 
obtained with these ne wer models are generally lower 
than those based on the iBruzual fc Chariot! (120031) mod- 
els f or galaxies with ages > 200 Myr (jReddv et all 
I2010T) . The relativ e contribution of the TP-AGB phase is 
still debated (e.g iMuzzin et al.l l200H: iKriek et alj|2010t 
Melb ourne et al.l 12012ft. a nd we note that adopting the 



Bruzual fc Chariot! (!2003[ ) models does not significantly 



alter our results. 

If such measurements were available, we corrected the 
broadband photometry (optical and/or AT s -band) for the 
effect of Lya emission/absorption and/or Ha emission. 
We did not explicitly correct for [OIII] emission, which 
lies in the AT s -band at z ~ 3 and the F160W band 
at z ~ 2, as [OIII] line measurements were not avail- 
able. However, neglecting the correction for Lya, Ha, 
and/or [OIII] emission for most of the galaxies in our 
sample results in differences in stellar masses and ages 
that are substantially smaller t han the marginali zed er- 
rors on ages and stellar masses (|Reddv et al.il2010n . This 
is due in part to the inclusion of the IRAC data where 
line contamination is not an issue, and where such data 
provide an additional lever arm to measure the strength 
of the Balmer and 4000 A breaks. 

For each galaxy, we considered a CSF model and ex- 
ponentially declining star formation histories with char- 
acteristic timescales r d = 10, 20, 50, 100, 200, 500, 1000, 
2000, and 5000 Myr. For comparison, we also investi- 
gated the effect on the stellar population parameters if 
we adopt exponentially rising star formation histories, 
where the SFR, ^P, is expressed as: 



(i) = * exp(i/T r ), 



(1) 



where ^ Q is the normalization factor, t is time (or age), 
and r r is the exponential timescale for the rising history. 
We have considered exponential timescales r r = 100, 200, 
500, 1000, 2000, and 5000 Myr. These star formation 
histories mimic linearly increasing ones if t agc <C T T . We 
further considered a range of ages spaced roughly log- 
arithmically between 50 and 5000 Myr, excluding ages 



G 



older than the age of the universe at the redshift of 
each galaxy. The lower age limit of 50Myr is adopted 
to reflect the dynamical timescale as inferred from ve- 
locity dispersion and size measurements of z ~ 2 LBGs 
(|Erb et alJl2006at lLaw et al.ll2009l 120121 ): the imposition 
of this age limit precludes galaxies from having unrealis- 
tic ages that are substantially younger than the dynam- 
ical timescale. As discussed below, we also investigate 
the effect of relaxing this age constraint and show how 
adopting a lower age limit (combined with a different at- 
tenuation curve) can resolve the discrepant SED-infcrrcd 
SFRs of young galaxies relative to those obtained from 
direct measurements of the SFRs derived from combin- 
ing UV and Spitzer/MlPS 24 /im data. F inally, redden- 
ing is taken into account by employing the lCalzetti et all 
(2000) attenuation curve (but see below) and allowing 
E(B — V) to range between 0.0 and 0.6. The choice 
of the Calzetti model is motivated by the good agree- 
ment between the Calzetti dust corrections and those de- 
termined from Spitzer/MlPS 24 /mi and Herschel/PACS 
100 and 160 p inferences of t he infrared luminosities 
(jReddv et al.ll2006bL [2010l I20B ) . 

The model SED at each r and age (i a gc) combination 
is reddened, redshifted, and attenuated blueward of rest- 
frame 1216 A for the opacity of the IGM using the lMadaul 
(1995) prescription. The best-fit normalization of this 
model is determined by minimizing its x 2 with respect to 
the observed U n GK+ JX S +F160W+1RAC (3.6-8.0 /mi) 
photometry. This normalization then determines the 
SFR and stellar mass. The model (and normalization) 
that gives the lowest x 2 is taken to be the best-fit SED. 
Typically there are several best-fit models that may ad- 
equately describe the observed photometry, even when 
the redshift is fixed to the spectroscopic value, though 
there is generally less variation in stellar mass than in the 
other para meters (r, t^, E(B — V)) among these best- 
fit models (iSawicki fc Yesd 119981 : iPapovich etaD 120011 : 
IShaplev et al.l l2001. 2003). Below, we consider a variety 
of star formation history models with different assump- 
tions for the age limit and attenuation curve (Table [3]). 

4. MULTI- WAVELENGTH CONSTRAINTS ON THE SFRS 
AND STELLAR POPULATIONS OF HIGH- REDSHIFT 
GALAXIES 

In this section, we compare the SFRs derived from SED 
fitting (SFR[SED]) with those calculated from combining 
the UV and MIPS data (SFR[IR+UV]). As discussed in 
Section[2j there are 302 galaxies in our sample with MIPS 
24 /mi observations and spectroscopic redshifts 1.5 < z < 
2.6 (121 of these galaxies are detected individually at 
24 /mi) ; it is at these redshifts where the 24 /mi fluxes 
are sensitive to the rest-frame 8 /mi emission, which in 
turn can be converted t o £tr . Using the procedure de- 
scribed in iReddv et al.l (|2010l ) , we fc-corrected the 24 /mi 
fluxes to estimate rest- frame 8 /mi luminosities (L s ). We 
then assumed a ratio Lir/Ls = 8.9 ± 1.3, as determined 
from a stacking analysis of the Herschel/PACS 100 and 
160 /mi data for a subset of the same galaxi es considered 
here (i.e., th ose in the GOODS -North field: IReddv et al.l 
120121 ). The iKennicuttl (|1998| ) relations are then used 
to convert Ltr and Ltjv (uncorrected for extinction) to 
SFRs, the sum of which gives the bolometric SFR. 

Note that bolometric SFR computed in this way is not 
entirely independent of the stellar population because the 



Table 3 

Description of SED Models 



Model Assumptions RMS a 



Model 


A 


Constant Star Formation 
Calzetti Attenuation Curve 
No Age Limit 


0.44 


Model 


B 


Declining Star Formation 
Calzetti Attenuation Curve 
No Age Limit 


0.48 


Model 


C 


Rising Star Formation 
Calzetti Attenuation Curve 
No Age Limit 


0.44 


Model 


D 


Constant Star Formation 
Calzetti Attenuation Curve 
iagc > 50 Myr 


0.44 


Model 


E 


Declining Star Formation 
Calzetti Attenuation Curve 
iagc > 50 Myr 


0.48 


Model 


F 


Rising Star Formation 
Calzetti Attenuation Curve 
iagc > 50 Myr 


0.46 


Model 


G 


Constant Star Formation 
Calzetti/SMC Attenuation Curves' 3 
iagc > 50 Myr 


0.46 


Model 


H 


Declining Star Formation 
Calzetti/SMC Attenuation Curves' 5 
iagc > 50 Myr 


0.50 


Model 


I 


Rising Star Formation 
Calzetti/SMC Attenuation Curves' 5 
iagc > 50 Myr 


0.48 



a RMS dispersion (in dex) about the best fit relation be- 
tween MIPS+UV SFR and SED inferred SFRs, as deter- 
mined from the expectation maximization (EM) paramet- 
ric estimator, for typical galaxies with Z/bol < 10 12 Lq and 

iagc > 100 Myr. 

b In this model, we have assumed the Calzetti et al. (2000) 
attenuation curve for galaxies with Calzetti-derived ages of 
iag a e' z > 100 Myr. For those galaxies with t^ z < 100 Myr, 
we remodeled their photometry assuming an SMC attenu- 
ation curve. 

conversion from L\jy to SFR will depend on the star for- 
mation history and age. This dependence is discussed in 
detail in Appendix For exponentially declining star 
formation histories where t > r^, the factor to con- 
vert Luv to SFR will be substantially smaller than the 
IKennicuttl (|1998l) value. More generally, for all of the 
models considered here (declining, constant, and rising), 
the conversion factor is at least a factor of two larger than 
the canonical value for ages < 10 Myr. For convenience, 
in Appendix [X] we provide formulae to compute the con- 
version factor between I/tjv and SFR for different star 
formation histories and ages, a ssuming the CB 11 stellar 
population synthesis models, a iSalpete rj ([1955J) IMF be- 
tween 0.1-100 M Q , and solar metallicity. Finally, we note 
that the factor to convert Lir to SFR is also some what 
dependent on galaxy age, with the IKennicuttl (|1998f) con- 
version valid for starbursts with ages < 100 Myr. At 
older ages, there are the competing effects of lower dust 
opacity and heating from older stellar populations that 
can modulate the conversion factor between Ltr and SFR 
(|KennicuttJ 1 19981 ) . These effects at older ages are likely 
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Figure 2. Comparison between SFRs derived from SED fitting and those computed from combining MIPS 24 /im and UV data. Points 
are color coded to differentiate galaxies with large bolometric luminosities (Lbol > 10 12 Lq), typical galaxies that form the bulk of the 
sample (Lbol < 10 12 Lq and t age > 100 Myr), and young galaxies with t agc < 100 Myr. Upper limits are points denote galaxies with 24 /im 
nondetections and detections, respectively. The large pentagons show the average 24 /tm plus UV determined SFR in bins of SED-inferred 
SFR for typical galaxies, as determined from a 24 fim stacking analysis (see the text). The results for a constant and rising star formation 
history, with no lower limit imposed on the age, are shown in the left and right panels, respectively. 

to be negligible for most of the high-redshift galaxies 
studied here due to their larger SFRs, and lower stellar 
masses, compared to the local galaxies used to calibrate 
the relationship between Lm and SFR. 

For simplicity, we assume the canonical Kcnnicutt 
(1998) relations to convert Lir and Ltjv to SFR, and 
discuss below how changing the Ltjv-SFR conversion fac- 
tor affects our results. The MIPS+UV derived SFRs are 
compared to the SED-inferred SFRs in Figure [5] Below, 
we discuss in turn the three sets of objects that differ- 
entiate themselves in the plane of SFR[IR+UV] versus 
SFR[SED]. 



SFR[IR+UV] and SFR[SED] of 0.44 dex (Table I3R) . 
Based on the stacking analysis and the survival analysis, 
we conclude that there is a good agreement between 
the MIPS+UV and SED derived SFRs for typical 
star-forming galaxies at z ~ 2. These results are not 
surprising because pr evious investigations have shown 
that on average the [Calzetti et all (|2000|) corrections 
for dust obscuration (which are assumed in the SED 
fitting procedure) reproduce the values estimated from 
mid and far-infrared, radio, and X-ray da t a for galaxies 
at z ~ 2 - 3 (e.g.. iNandra et all l200l ISeibert et all 
2001 iReddv fc Steidell 120041: IReddv et all 12006 



Daddi et all 120071: IPannella et all 120091: IReddv et all 
20101: IMaedis et alJl2o!onReddv et al.H2012D . 



4.1. SFRs of Typical Star-forming Galaxies at z ~ 2 

This subsection highlights the results found for typical 
star-forming galaxies, defined as those galaxies with best- 
fit ages > 100 Myr and bolometric luminosities L^ \ < 
10 12 Lq. Here, -L do i 1S taken as the sum of the UV and IR 
luminosities, where Lir is computed from the 24 /im data 
using the Ltr/Ls ratio discussed above. Because 60% of 
the galaxies are not detected at 24 /im, we have stacked 
the 24 /im data for typical galaxies in bins of SFR[SED] 
to determine the average bolometric SFR. The stacking 
procedure is discussed in IReddv et all (|2006bD . Uncer- 
tainties in the stacked bolometric SFR are computed by 
combining in quadrature the measurement uncertainty 
in the mean 24 /im flux and the uncertainty in the mean 
UV luminosity of galaxies contributing to the stack. 

SFR[IR+UV] derived in this manner agrees well on 
average with SFR[SED] assuming constant star forma- 
tion and no age limit (large pentagons in the left panel 
of Figure [5]) . Taking into acco unt upper limits u sing the 
ASURV statistical package (jlsobe et al.l I1986D . which 
includes the expectation maximization (EM) parametric 
survival estimator for censored data, we compute an 
rms dispersion about the best-fit linear relation between 



4.2. SFRs of Ultraluminous Infrared Galaxies at z ~ 2 

Turning our attention to ultraluminous infrared galax- 
ies (ULIRGs), we find that those galaxies with Lhoi > 
10 12 Lq have bolometric SFRs that exceed by up to a 
factor of ten the SFRs [SED]. The 24 /im fluxes of the 
IR-luminous sources in our sample tend to over-predict 
their L\n by a factor of w 2, relative to the IR estimates 
obtained by includ ing far-IR data (e.g., from Herschel; 
IReddv et all I2012I). Further, as shown in several other 
investigations, the lCalzetti et all (|2000D dust corrections 
for such objects are typically too low due to the fact that 
much of the star formation is completely obscured by 
dust, and hence the UV slope decouples from extinction 
for such highly obscured galax ies (e.g., iGoldader et all 
l2lMlReddv et al.ll2006bl l2010h . 

4.3. SFRs of Young Galaxies at z ~ 2 

A noted problem in stellar population modeling is the 
distribution of unrealistically young ages derived for non- 

13 The RMS values listed in Table [3] merely indicate the rms 
about the best-fit relation between SFR[SED] and SFR[IR+UV], 



and arc not meant to indicate the 
two quantities. 



goodness of fit" between the 



negligible fractions of galaxies in high-redshift samples, 
particularly tho se selected in the rest-UV or rest-optical 
(Sh aplev et aHl200ll: IMaraston et al.|[20loh . This is com- 
monly referred to as the "outshining" problem, where the 
SED of a galaxy may be dominated by the youngest stel- 
lar population even at near-IR wavelengths, in which case 
the best-fit models (irrespective of the star formation his- 
tory) may have very young (< 100 Myr) ages with very 
large values of SFR and reddening (e.g.. IMaraston et al.l 
120101: iWuvts et all 1201 If ). We have investigated the va- 
lidity of these young models by cross-checking the SFRs 
derived using them with those measured directly from 
MIPS 24 fim and UV data. 

For such young (t age < 100 Myr) objects (the third set 
of objects denoted in Figure |2J), we find that SFRs[SED] 
are systematically larger than SFRs[IR+UV]. This off- 
set cannot be completely accounted for by a change in 
the conversion between UV luminosity and SFR. In par- 
ticular, the bulk of these young galaxies have ages (de- 
rived with the CSF assumption) as short as 7 Myr. For 
these ages, the conversion between UV luminosity and 
SFR is about a factor of two larger than in the CSF 
case where t agc > 100 Myr (Appendix Given the 
typical dust attenua tion of these "young" gal axies of 
Lm/Lvw w 1 (e.g.. iReddv et al.l [2006H 120121) . modi- 
fying the UV-SFR conversion upward by a factor of two 
will increase SFR[IR+UV] by 50%, or 0.18 dex. This dif- 
ference is not sufficient to account for the offset between 
SFR[SED] and SFR[IR+UV] for the young subsample. 

Galaxies at z ~ 2 that are identified as being young 
based on their CSF model fits have UV slopes, (3, which 
are on average redder than the slopes of older galaxies at 
the same redshifts, with a difference in /3 of (8/3) ~ 0.3 
(Figure [3]). If we assume that a Calzetti attenuation 
curve applies for such galaxies, then we would find them 
to be more dust obscured, given their redder /?, and 
hence to have larger bolometric SFRs, than the "typ- 
ical" galaxies discussed above. Clearly, the very large 
dust obscuration and bolometric SFRs inferred for such 
galaxies cannot be correct, based on the comparison of 
these SFRs with those derived directly from combining 
the MIPS and UV data (Figure EJ. 

4.4. Results for Other Star Formation Histories 

We conclude this discussion by noting that the assump- 
tion of exponentially rising star formation histories re- 
sults in the same qualitative behavior for the three sets 
of galaxies discussed above (typical galaxies, ULIRGs, 
and young galaxies; Figure [5]). The RMS around the 
best-fit relation between SFR[SED] and SFR[IR+UV] is 
listed in Table [3] In contrast, the assumption of ex- 
ponentially declining star formation histories results in 
SFRs [SED] for typical galaxies that lie systematically 
below SFRs[IR+UV] (Figure |4j. A similar conclusion 
is reached when comparing dust-corrected Ha estimates 
of the SFR with those computed from fitting exponen- 
tially declining models to the broa dband photometr y of 
a subset of galaxies in our sample (|Erb et al.ll2006aj ). 

The discrepancies that arise from adopting an expo- 
nentially declining star formation history stem from two 
effects. The first is that the redness of the UV contin- 
uum is, to a greater extent, attributed to older stars in 
the case when t^/r^ > 1, relative to the CSF case. 
Hence, the E(B — V) for declining models will on aver- 




Figure 3. /3 distribution for galaxies with Calzetti-derived ages 
> 100 Myr (red histogram) compared to that of galaxies with 
Calzetti-derived ages < 100 M yr (blue histogram) , in the sample 
of 392 UV selected galaxies offRcddy et al. (201(jj). The ages as- 
sume a CSF model. A two-sided K-S test indicates a probability 
P < 0.01 that the two distributions are drawn from the same par- 
ent sample. On average, the younger galaxies have UV slopes that 
are (80) ~ 0.3 redder than those of older galaxies. 

age be lower, implying lower SFRs. The second effect 
is that the ratio of the SFR to UV continuum is gener- 
ally lower for galaxies w ith smaller Td (Appendix [X] see 
also lShaplev et afl l2005 for a discussion of these points), 
which also results in sma l ler SF Rs for a given UV lu- 
minosity. iShaplev et al.l (j2005[ ) discuss other reasons 
why extreme declining star formation histories where 
r c i < 100 Myr and tage/^d 3> 1 are implausible for the 
vast majority of galaxies in our sample. In particular, the 
correspondence be tween the Ha and UV derived SFR s 
for z ~ 2 galaxies (|Erb et al.ll2006at IReddv et alJIMot l. 
as well as the presence of O star features in the UV 
spectra of galaxies in our sample — e.g., Lya emission, 
and Wolf Rayet signatures including broad Hell 1640 A 
and P-Cygni wind features in the CIV and SilV lines 
(jShaplev et alj2003t iQuider et al.ll2009fl — imply that the 
UV continuum is likely dominated by O stars, contrary 
to the expectation if i age /Td ^> 1. 

The importance of the comparison to the MIPS+UV 
SFRs is underscored by the fact that one cannot dis- 
criminate between these various star formation histories, 
even the simple ones considered here, from the broad- 
band SEDs alone. As the right panel of Figure @] shows, 
the x 2 distributions of the best-fit SEDs for the different 
star formation histories are roughly similar. The MIPS 
data allow us to break this impasse, and they can be used 
to demonstrate clearly that, on average, the exponen- 
tially declining models yield SFRs that are statistically 
inconsistent with those obtained from direct measures of 
the SFRsEB Given this finding, we focus the subsequent 
discussion on constant and rising star formation histo- 

14 Wuyts et al. (2011) show that restricting declining star for- 
mation histories to e-folding times > 300 Myr result in a better 
agreement between SED-inferred SFRs and those obtained from 
combining UV and IR data. Doing the same for our sample, we 
obtain a median value of t/r ~ 0.4 which implies a behavior closer 
to that of the CSF models. 
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Figure 4. (Left:) Same as Figure[2] where we have assumed the best-fit exponentially declining star formation history. (Right:) Histograms 
of the x 2 distributions of the best-fit SED models to the broadband photometry for 1959 galaxies with spectroscopic redshifts 1.4 < z < 3.7, 
assuming declining, constant, and rising star formation histories (corresponding to Models A, B, and C, respectively, as listed in Tabled- 
The histograms have been normalized to the total number of galaxies (1959). 
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Figure 5. Ratio of SED-inferred SFR and SFR derived from com- 
bining MIPS and UV data, vs. stellar population age, for Model C 
(Table [3] rising star formation history with no age limit). Symbols 
are the same as in Figure [2] The large red pentagon denotes the 
median value based on stacking the 24 fim data for the galaxies 
with t ago < 100 Myr. For Model C, the SED-inferred SFRs of the 
young galaxies exceed by an order of magnitude those estimates 
from combining the MIPS and UV data. 



RESOLVING THE CONFLICTING SFRS FOR YOUNG 
GALAXIES 



As noted in Section FOl all of the star formation histo- 
ries considered here predict SFRs that are substantially 
larger than those computed from the MIPS and UV data 
for young galaxies with i agc < 100 Myr. Figure [5] shows 
that this systematic difference is a strong function of 
age, being most severe for galaxies with t age < 50 Myr. 
Given that the true ages are unlikely to be significantly 
smaller than the dynamical timescale of ~ 50 Myr, we 
have imposed the restriction £ ago > t^yn — 50 Myr when 
fitting the SEDs. Note that the exact limit in age 



is unimportant given the relatively few galaxies with 
20 < £ ag c ^ 100 Myr. Adopting an age restriction, the 
comparison between SED and MIPS+UV SFRs is shown 
in Figure [6] At face value, based on the individual data 
points, the agreement between the SFRs for galaxies with 
iago < 100 Myr is better once we have restricted the age 
to be older than £ ago = 50 Myr. However, a MIPS stack 
for those galaxies with 50 < t agc < 100 Myr results in 
a formal non-detection at 24 /im, inconsistent at the 3cr 
level with our expectation based on the SFRs [SED] for 
these galaxies (red pentagons in Figure [6|). 

5.1. Age Dependence of the UV Attenuation Curve 

One possible solution to this discrepancy comes from 
the correlation between the dust obscuration, Lir/LtjVj 
and UV slop e (0) for these young galaxies. In particular, 
iReddv et abl (|2006bl l2010h found that such galaxies have 
redder UV slopes for a given dust attenuation relative to 
older galaxies with i agc > 100 Myr. This effect may stem 
from a geometrically different distribution of dust with 
respect to the stars in a galaxy where, for young galaxies, 
this dust may preferentially lie foreground to the stars. 
Regardless, the attenuation curve obtained for the young 
galaxies (on average) appears to mimic that of the Small 
Magellanic Cloud (SMC) attenuation curve. A steeper 
attenuation curve (e.g., SMC curve) will result in older 
ages than those obtained in the Calzetti case; because 
not as much dust is needed to redden the continuum for 
a foreground screen of dust relative to a uniform distri- 
bution of dust and stars, the redness of the UV slope 
will be attributed more to an older stellar population 
than dust attenuation. If we remodel those galaxies with 
Calzetti-inferred ages 50 < i agc < 100 Myr with the SMC 
attenuation curve, then the fraction of galaxies that are 
still considered to be "young" (£ agc < 100 Myr) is reduced 
by at least a factor of nine. In the CSF case, the number 
of galaxies with t age < 100 Myr drops from 36 to just 
4. In the exponentially rising case, this number drops 
from 36 to just 2. The small fraction of galaxies that 
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Figure 6. Same as Figure [2] where we have imposed a lower limit on the age of t agc = 50Myr in the SED modeling including constant 
(left) and rising (right) star formation histories. The large red pentagon denotes the stacked results for galaxies with 50 < t age < 100 Myr. 
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50 Myr in the SED modeling, and we have remodeled 



rising (right) star formation histories. The green arrow in the left panel indicates the 3 it upper limit to SFR[IR+UV] for the 4 24 /im 
undetected galaxies that are considered young under the assumption of the SMC attenuation curve. 

At first glance, the adoption of a lower limit in age 
that is equivalent to the dynamical timescale, combined 
with a different attenuation curve for the "young" galax- 
ies, may appear to be a contrived and inelegant solution 
to resolving the discrepancy between SFR[IR+UV] and 
SFR[SED] for such galaxies. However, our primary goal 
is to derive SED parameters that are based on phys- 
ically motivated ages (i.e., dynamical time constraints) 
and are consistent with obser vations of the dustiness an d 
UV slopes for such galaxies (jReddv et al.ll2006bl 



are still considered young under the SMC assumption 
have upper limits in SFR[IR+UV] that are consistent 
with SFR[SED] when we assume an SMC attenuation 
curve (Figure EJ). Hence, if we assume an SMC attenua- 
tion curve, there is no longer a conflict between the upper 
limit to SFR[IR+UV] (from the MlPS-nondetection) and 
SFR[SED] for the young galaxiesE] 



15 As we discuss in Section l6.2l there are certain star formation 
histories (namely rising ones) where the age of galaxy becomes 
an ill-defined quantity, thus obfuscating the distinction between 
"young" and "old" galaxies as discussed here. In these cases, it 
is useful to distinguish galaxies based on their stellar masses (or 
metallicities, if such measurements are available). From the SED- 
fitting assuming a CSF history, the vast majority of those galaxies 
with a ges < 100 Myr also have M* < 10 Mq. As wc show in 
Section [62] the assumption of a rising star formation history does 
not significantly alter the stellar masses relative to those obtained 



12010) . 

A consistent treatment of the multi-wavelength data al- 
lows us to resolve the discrepancy between the different 



under the CSF case and, as such, one can just as easily adopt 
the stellar mass threshold of M* 10 10 Mq to reach the same 
conclusions regarding the validity of the various dust attenuation 
curves for z ~ 2 galaxies. 
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measures of SFRs for young galaxies. 

5.2. Summary of SFR Comparisons 

In the previous sections, we presented the comparison 
of the SFRs determined from independent indicators of 
dust (MIPS 24 /im) to those computed from fitting the 
broadband SEDs of galaxies at redshifts 1.5 < z < 2.6. 
Three sets of galaxies are identified in this comparison: 
(1) ULIRGs with L hoi « L m > 10 12 L Q ; (2) typical star- 
forming galaxies with Lboi < 10 12 L Q and ages t agc > 
100 Myr; and (3) "young" galaxies with L ho i < 10 12 L Q 
and ages t agc < 100 Myr. For typical galaxies, an ex- 
ponentially declining star formation history yields SFRs 
that are inconsistent with those obtained from combin- 
ing the MIPS and UV data. Assuming constant or rising 
star formation histories yields SFRs that are in reason- 
able agreement with the MIPS+UV determinations. Al- 
ternatively, none of the simple star formation histories 
considered here (constant, declining, or rising) are able 
to reproduce the lower SFRs found for "young" galax- 
ies. We explore several possibilities for resolving this 
discrepancy, including dynamical time constraints and 
systematics in the UV attenuation curve. A physically 
plausible solution — and one which is consistent with 
independent measurements of the a ge dependence of the 
dust attenuation curve at z ~ 2 — 3 (Red dv et al.l l2006b. 
12010ft — is to adopt a lower limit in age corresponding to 
the dynamical timescale typical of galaxies in our sam- 
ple (~ 50 Myr), and to assume a steeper (e.g., "SMC") 
attenuation curve for the young galaxies (corresponding 
to Models G and I in Table [3]). Based on these findings, 
we proceed in the next section with a discussion of the 
stellar masses and ages derived from constant and rising 
star formation histories. 

6. INTERDEPENDENCES OF THE SED PARAMETERS 

Having just discussed the most plausible star forma- 
tion histories that are constrainable with the available 
data, in this section we focus on the systematic varia- 
tions in age and stellar mass with the overall form of 
the star formation history (constant versus rising), the 
typical uncertainties in ages and stellar masses, and the 
correlation between SFR and stellar mass. 

6.1. Systematic Variations in Age and Stellar Mass 

The ages and stellar masses derived for 1959 objects 
with spectroscopic redshifts 1.4 < z < 3.7, assuming con- 
stant and exponentially rising star formation histories, 
are shown in Figure [5] Histograms of the age and stellar 
mass distributions are presented in Figure [9] The stel- 
lar mass for any given galaxy derived with a rising star 
formation history is essentially identical within the un- 
certainties (that stem primarily from photometric errors 
and the degeneracy between the parameters being fit) to 
those obtained with a CSF history. The primary differ- 
ence, therefore, is in the distribution of best-fit ages; ages 
derived with a rising star formation history are older by 
w .12dex. The fraction of galaxies in our spectroscopic 
sample with derived ages older than 1 Gyr increases from 
20% to 32% when we assume rising star formation histo- 
ries versus CSF. 

We can better understand these trends by examining 
the evolution in the stellar mass-to-light ratio (M/L) as a 



function of galaxy age for the different star formation his- 
tories considered here. At early times, for ages i agc < r r , 
the light is dominated by current star formation and the 
M/L ratio is low. As the stellar mass builds up, older 
stars will begin to contribute significantly to the light, 
causing the M/L ratio to increase. Neglecting dust red- 
dening, for t agc 3> r r , the M/L ratio stabilizes to some 
constant value. In reality, the M/L ratio will still in- 
crease (even though the ratio of stellar mass to total 
SFR may remain roughly constant), because dust red- 
dening will increase with SFR, and hence the bolomet- 
ric light will suffer an increasing amount of obscuration 
(Figure HDJ). 

Regardless, as shown in Figure [TOl the M/L ratio (at 
rest- frame 1.1 /im, corresponding roughly to an observed 
wavelength at z ~ 2.3 that coincides with IRAC chan- 
nel 1) increases more slowly with age for galaxies with 
smaller values of r r . Hence, for the same M/L ratio, ris- 
ing star formation histories correspond to an older age 
than a CSF history. Note that these trends hold on av- 
erage, and there is still a fair fraction of galaxies where 
the ages (and stellar masses) do not change substantially 
between the CSF and rising star formation histories. In 
particular, approximately 55% of galaxies in our sample 
have a best-fit r r = 5000 Myr, and it is for these longer 
exponential timescales that the M/L ratio is not substan- 
tially different than that obtained in the CSF case when 
iagc < t t . In any case, the stellar mass distributions do 
not change substantially when we adopt a rising versus 
CSF history, and in both cases, the mean stellar mass of 
galaxies in our spectroscopically confirmed (1Z < 25.5) 
sample at z ~ 2.3 is (M») = 1.6 x 10 10 M Q . As we dis- 
cuss in Section 17.2. 1[ the mean stellar mass is a strong 
function of absolute magnitude, and including galaxies 
fainter than our spectroscopic limit will result in a lower 
mean stellar mass. 

6.2. The "Ages" of Galaxies with Rising Star 
Formation Histories 

As just discussed, rising star formation histories re- 
quire older ages than CSF histories to achieve a given 
mass-to-light ratio. In general, the ages derived with ei- 
ther a constant or rising star formation history will be 
a lower limit to the true age, because of the possibility 
that there may be an underlying older stellar population 
whose near-IR light is overwhelmed by the near-IR light 
from current star formation (i.e., the "outshining" prob- 
lem as discussed above) . A further systematic effect per- 
tains to "ages" derived under the assumption of different 
star formation histories. In the simplest case for constant 
star formation, the age is simply determined by the time 
required to build up the observed stellar mass given the 
current rate of star formation. Integrating Equation Q] 
yields the stellar mass after time t for an exponentially 
rising star formation history: 

M.(t) = * T r [exp(t/r r ) - 1], (2) 

For simplicity, we have ignored the gas recycling fraction 
(i.e., the fraction of gas released back into the interstellar 
medium (ISM) from supernovae and stellar winds), as 
adopting this correction will simply adjust downwards 
the st ellar mass by some multi plicative factor close to 
w 0.7 (|Bruzual fc Chariot! [2003h . With this parameteri- 
zation, the stellar mass M* = M at time t = 0. In this 
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Figure 8. Comparison of the ages (left) and stellar masses (right) between the constant (Model G) and rising star formation models 
(Model I), respectively. 

the currently observed value with a rising and CSF mod- 
els. Of course, for a given currently observed SFR, the 
rising history (as compared to a CSF model) will implies 
a longer amount of time to build up the current stellar 
mass given that the SFR was lower in the past. In Sec- 
tion [51 we discuss the mass doubling time, derived ages, 
and implied formation redshifts, in the context of direct 
observations of z > 3 galaxies. 
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Figure 9. Histograms of age and stellar mass for the constant 
(red) and rising (blue) star formation histories. The vertical lines 
in the left panel indicate the mean age found when assuming the 
constant and rising star formation histories. 

case, the "age" is well defined in a mathematical sense, 
and these are the ages that we have referred to in the 
previous sections. Note that the specific SFR 



m = * (t) 



exp(t/r r ) 



M*(i) Tr[exp(t/T r )-1] 



(3) 



6.3. Uncertainties in Age and Stellar Mass with Simple 
Star Formation Histories 

As just discussed, the assumption of rising versus CSF 
histories results in systematic differences in the stel- 
lar population ages. Other systematics affecting ages 
and stellar masses include uncertainties in the IMF, 
differenc es between stellar popu l ation synthesis mod- 
els (e.g.. iBruzual fc Charlotll2003l; IMaraston et alj|2006t 
CB11), and the degree to which simple star forma- 
tion histories can capture the complexi ty of the "real" 
star f ormation history of a galaxy (e.g.. [Papovich et al.l 



is roughly constant as a function of time if exp(t/r r ) 3> 1. 
As long as this condition is satisfied, the age t can be 
arbitrarily old. In that case, we would simply adjust 
the normalization \l/ for a given t in order to recover 
the observed stellar mass. As a simple example, if we 
assume that r r » 400 Myr and we consider ages of t = 1 
and 2Gyr, then we must set « 3.6 and 0.3M Q yr _1 , 
respectively, in order for the current stellar mass to be 
M* w 1.6 x 10 10 Mq. Given the ambiguity of the derived 
"age" for a rising star formation history, it is useful to 
define a characteristic timescale t c over which the galaxy 
has doubled its stellar mass: 



20011: iShapley et al ll200l 120051: iMarchesini et al.l 1200 
Muzzin et al.ll2009t IMaraston et al.ll2010UPapovich. et alJ 



tr. 



1 



1/2 



r r ln(2) = - ln(2) 



(4) 



where t±/2 indicates the time at which the galaxy had 
half of its current stellar mass. We show the compar- 
ison between t c and t age /2 for the rising and constant 
star formation models, respectively, in Figure [TT] The 
comparison shows that there is a better correspondence 
between the times required to double the stellar mass to 



20111) . While the uncertainties in the complexity of the 
star formation history of any individual galaxy may be 
large, if such variations are stochastic (e.g., such as hav- 
ing multiple independent bursts of star formation), then 
large galaxy samples can be used to average over these 
stochastic effects and yield insight into the character- 
istic, or typical, star formation history. In addition 
to these systematic effects, there are also uncertainties 
related to our ability to accurately constrain the stel- 
lar population model given the flux errors and limited 
wavelength coverage of photometry. These effects result 
in typical fractional uncertainties of er agc /(£ agc ) = 0.5 
and ayi*/{M*) = 0.4, re spectively, as determined from 
Monte Carlo simulations (|Shaplev et al.| [2005: Er b et all 
l2006bD . 

6.4. Relationship between Star Formation Rate and 
Stellar Mass 

The relationship between SFR[IR+UV] and stellar 
mass for the 302 galaxies at 1.5 < z spcc < 2.6 with MIPS 
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Figure 10. Left: Stellar mass to light ratio at 1.1 /an (in solar units) as a function of age for CSF and exponentially rising star formation 
histories, with r r = 100 to 5000 Myr. The curves assume E(B — V) = 0. Right: Evolution in M/L ratio including the effects of dust 
attenuation (dashed lines), where we have assumed the relationship between SFR and E(B — V) as determined from the SED modeling, 
compared to the cases with no dust attenuation, E(B — V) = (solid lines). For this example, we consider r r = 1000 and 2000 Myr, and 
an initial SFR of 20 Mq yr" 1 . The difference between the curves for r r = 1000 Myr is larger than the difference between the T r = 2000 Myr 
curves because the SFR, and hence dust reddening, rises more rapidly for smaller values of r r . 

of « 0.37 dex, and a slope of 0.97 ± 0.05. 

Because the relationship between SFR and stellar mass 
has been s tudied extensively over a large range of red- 
shifts (e.g..lBrinchmann etaT1l2004t IReddv et al]|2006b t 
Noeske et all 120071 iDaddi et all 120071 : iStark et al.ll200£ 




Gonzalez et al . 2010; iSawickj II2011D . it is useful to ex- 



15/2 (Myr) 

Figure 11. Comparison between the characteristic timescalcs re- 
quired to double the stellar mass to the present value with a rising 
(t c tslng ) and constant (t^, F '/2) star formation histories for the 
full SED sample of 1959 galaxies. 

24 /im observations (121 of which are detected at 24 /im) 
is shown in the left panel of Figure Q21 Rising star forma- 
tion histories are assumed; stellar masses derived from a 
CSF history will be similar. As found by several previ- 
ous investigations at similar redshifts (e.g.. IReddv et all 
I2006bt iPapovich et all [20061: IDaddi et all 120071 ). we find 
a relatively tight positive correlation between SFR and 
stellar mass (though we discuss below potential biases), 
which has typically been interpreted as an indication that 
most of the stellar mass at these redshifts accumulated 
in a relatively "smooth" manner and not in major merg- 
ers. Taking into account the upper limits with a survival 
analysis, we find an rms scatter about the best-fit linear 
relation between the log of SFR and log of stellar mass 



amine how sample selection may bias the measurement 
of this relationship. The right panel of Figure [12] shows 
the correlation between SFR and M*, where SFR is de- 
termined solely from the SED fitting (using Model I), 
for a larger sample of 1244 galaxies with spectroscopic 
redshifts in the same range as the MIPS sample. For 
comparison, the distribution of SFR versus stellar mass 
for the subset of 705 galaxies with at least two detec- 
tions longward of the Balmer and 4000 A breaks is also 
shown. We note that the distribution for this subset is 
similar to that of the larger sample, implying that the 
inclusion of galaxies with just a single point longward of 
the Balmer and 4000 A breaks does not systematically 
bias the distribution of SFR versus M*. In the SED 
fitting procedure, both the SFR and stellar mass are de- 
termined by the normalization of the model SED to the 
observed photometry (Section [3|) . As such, SFR[SED] 
and .M* are highly correlated and exhibit an intrinsic 
slope of unity. However, a linear least squares fit to the 
data imply a much shallower slope of 0.30. This bias is 
a result of the fact that our sample is selected based on 
UV luminosity (e.g., and not stellar mass). Hence, there 
will be a Malmquist bias of selecting galaxies with larger 
SFRs at a given stellar mass for the lowest stellar mass 
galaxies in our sample. This bias is quantified via simu- 
lations that are presented in Appendix [Bj The effect of 
this bias on the best-fit relation between stellar mass and 
near-IR magnitude is discussed in Appendix [C] Finally, 
Appendix [E] examines the effect of these biases on the 
determination of mean stellar mass from a UV selected 
sample. 
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Figure 12. Left: Relationship between SFR, determined from combining MIPS and UV data, and stellar mass as determined from SED 
fitting for 302 spectroscopically confirmed galaxies at 1.5 < z < 2.6. We have assumed rising star formation histories to derive the stellar 
masses (e.g., using Model I [Table l3l ) . Upper limits and detections are indicated by the arrows and hollow circles, respectively, and the 
darker filled symbols and upper limit denote the stacked SFR in bins of stellar mass. The best-fit linear correlation between SFR and M* 
for the spectroscopic sample, taking into account upper limits, is shown by th e solid blue line. Th e large open stars are for photometrically 
selected faint z ~ 2 galaxies. The red line indicates the relationship found by Daddi ct al. (2001]) at similar redshifts. Right: Same as left 
panel, where SFRs are determined from SED fitting of a larger sample of 1264 galaxies with spectroscopic redshifts in the same range as 
the MIPS sample (cyan points). The purple points show the distribution for the subset of 705 galaxies that have at least two detections 
longward of the Balmer/4000 A breaks. The solid line shows the intrinsic correlation between SFR and M„, and the dashed line indicates 
a linear least squares fit to the data. 



A comparison of the left and right panels of Figure [12] 
shows that the Malmquist bias is much less noticeable 
when the SFR is determined independently of the stel- 
lar mass (which is not the case when fitting the SEDs 
with stellar population models), primarily because the 
SFRs[IR+UV] have additional scatter associated with 
them (related to the dispersion in dust extinction at a 
given UV luminosity) that is decoupled from the method 
we used to infer M* from the SED fitting. As we dis- 
cuss below, the near-unity slope of the SFR-M* relation, 
where SFRs are determined from the MIPS+UV data, 
appears to hold to lower stellar masses, and it implies 
that the specific SFR depends only weakly on stellar 
mass. We will return to the implication of this result 
in Section [H 

Before proceeding, wc comment briefly on a poten- 
tial bias that may exist at the bright, high-mass end 
of the SFR-M* re lation. As noted in Section SI the 
IReddv et all ((20101) conversion between rest-frame 8 /im 
and total IR luminosity reproduces the average Ltr, com- 
puted using st acked Herschel/PACS 100 and 160 /im data 
(jReddv et all 120121 ). This conversion is also shown to 
over-predict Lir for the most bolometrically luminous 
galaxies in our sample. Correcting for this effect will 
lower SFR[IR+UV] by up to a factor of two, resulting in 
a shallower slope of SFR-M* at these high stellar masses 
where M* > 10 n MQ. On the other hand, at a given 
high stellar mass, our UV selected sample will be biased 
against the dustiest and hence most heavily star-forming 
galaxies. This is due the fact that such dusty galaxies 
either (a) have UV colors that are too red to satisfy the 
color criteria or (b) are too f aint to be represented i n 
the spectroscopic sample (e.g. . IReddv et a l. 2008, 2010). 
Inclusion of these missing dusty galaxies would have the 
opposite effect of steepening the slope of the SFR-M* re- 
lation at high stellar masses. Regardless of these two 
competing biases at the high stellar mass end, correcting 
for them is unlikely to alter significantly the overall slope 
of the SFR-M* relation, given that such high SFR galax- 
ies constitute only a small fraction of the number den- 



sity of L* and fainter galaxies (e.g.. iReddy et al.l 120081 : 
IReddv fc SteTdel[20CJ9tMagnelli et al.ll201lD . 

6.5. Summary of Age and Stellar Mass Distributions, 
and SFR-M* Relation 

Modeling the broadband photometry of spectroscopi- 
cally confirmed galaxies in the z ~ 2 — 3 sample allows 
us to examine the distribution of ages and stellar masses. 
We have shown that rising star formation histories gener- 
ally yield older ages for a given mass to light ratio relative 
to the ages obtained with a CSF history. More gener- 
ally, the "age" of galaxy is not well-defined for a rising 
star formation history, in the sense that the age can be 
varied simultaneously with the normalization of the star 
formation history to yield the same value of current SFR 
and stellar mass, as long as exp(tr) 3> 1. Using SFRs 
that are derived independently of the stellar masses, we 
find an intrinsic SFR versus stellar mass correlation with 
roughly unity slope at z ~ 2. So far, we have concerned 
ourselves with the SED sample which, by construction, 
only includes those galaxies that had at least one detec- 
tion longward of the 4000 A break. However, there exist 
substantial numbers of UV-selected galaxies that are not 
represented in the SED sample because they are unde- 
tected longward of the break. More generally, we expect 
the SED sample to be biased to galaxies with larger stel- 
lar masses at a given SFR relative to the UV sample as 
a whole. In the next section, we explore this bias by 
stacking the longer wavelength data for galaxies in our 
sample, and we also extend the results on the SFR-M* 
relation to UV faint galaxies at the same redshifts. 

7. MASS-TO-LIGHT RATIOS OF Z ~ 2 - 3 GALAXIES 

7.1. Rest-frame Near-IR M/L Ratios 

Because the data required to model the stellar popu- 
lations do not exist for every galaxy in our sample (e.g., 
those galaxies that did not have near-IR or IRAC cov- 
erage, as well as those undetected at these longer wave- 
lengths), it is useful to establish some empirical relation- 
ship between the luminosity of a galaxy and its stellar 
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Table 4 

Mass to Light Ratios 



Rest-frame A 
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Spearman's a h 


RMS C 




( M\ o 
^ "~F" ^ max 


0.50 


98 


5.5 


0.33 


0.07 


4.4 


0.67 


491 


14.5 


0.32 


0.06 


4.4 


1.1 


643 


22.3 


0.20 


0.14 


4.4 


1.4 


673 


23.0 


0.18 


0.18 


4.4 


1.8 


180 


11.8 


0.19 


0.24 


2.4 


2.4 


190 


10.9 


0.23 


0.21 


4.4 



a Number of spectroscopically confirmed galaxies. 

b Significance of correlation between log stellar mass and absolute mag- 
nitude computed from Spearman's rank correlation test, in units of stan- 
dard deviation. 

c RMS of data about best-fit linear relation between log stellar mass and 
absolute magnitude. 

d Minimum observed mass to light ratio in our sample, in units of the 
Sun. 

c Maximum observed mass to light ratio in our sample, in units of the 
Sun. 



mass. Given the stellar mass to light ratio (M/L), we 
can then estimate the stellar mass of galaxy based simply 
on the flux at a given wavelength. The relation between 
stellar mass and rest-frame near-IR light as probed by 
the IRAC observations is shown in Figure [T31 We have 
also computed these relations at rest-frame optical wave- 
lengths using the F160W and if s -band data As summa- 
rized in Table |4j the dispersions between the rest-frame 
near-IR luminosities and stellar masses are lower than 
those found between rest-frame optical luminosity and 
stellar mass, but we note the large variation in mass- 
to-light ratio even at rest-frame wavelengths where the 
stellar emission peaks, around 1.6 /im. The current star 
formation in a galaxy may outshine the light from the 
older stars at rest-frame optical and near-IR wavelengths. 
In these cases, the ages (and masses) typically reflect 
those of the current star formation episode, though two 
component models can be used to determine an upper 
limit to the hidden stellar mass in such galaxies (e.g., 
IShaplev et al.ll2005h . For the subsequent discussion, we 
focus on results using the IRAC 3.6 /mi data, noting that 
our results (e.g., derived stellar masses) do not change 
substantially when we use data from the other IRAC 
channels. 

Note that there is a systematic trend towards lower 
M/L at lower stellar mass or fainter near-IR luminosity; 
i.e., the best-fit relation between near-IR luminosity and 
stellar mass does not fall onto a line of constant M/L. 
At face value, this systematic trend implies that galaxies 
with lower stellar masses have larger specific SFRs (</>), 
resulting in lower M/L ratios. However, as discussed in 
Section 16.41 and Appendix [Bj there is Malmquist bias of 
selecting galaxies with larger SFRs, and hence larger cf>, 
for the lowest mass galaxies in our spectroscopic sam- 
ple. This directly affects the best-fit relation between 
stellar mass and near-IR luminosity, as discussed in Ap- 
pendix [Cj The basic conclusion from these different best- 
fit relations is that there are sufficient biases induced by 
Malmquist effects and not correcting for dust extinction 
that can produce an erroneously strong trend between 
M/L ratio and luminosity. For the subsequent discus- 
sion, we employ the relations discussed in Appendix [C] 
to convert the average near-IR magnitude to stellar mass 
for galaxies of a given UV luminosity. 
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Figure 13. Relationship between stellar mass and absolute mag- 
nitude measured in the four IRAC channels. Approximate rest- 
frame wavelengths are indicated in each panel, along with the num- 
ber of objects used to define the relationship. Short dashed lines 
indicate the minimum and maximum M/L ratio. 

7.2. Rest-frame UV Mass-to-light Ratio 

The scatter in the M/L ratio increases towards shorter 
wavelengths, as the current star formation dominates the 
emission. This can already be seen in the larger scatter in 
M/L at optical wavelengths relative to that found in the 
near-IR, and we would expect the maximum dispersion 
to occur at UV wavelengths. Here, we investigate the 
mean and dispersion in M/L ratio at UV wavelengths, 
with the aim of quantifying the stellar masses and star 
formation histories of UV faint galaxies. 

7.2.1. Trend Between Stellar Mass and UV Luminosity in 
the Spectroscopic Sample 

Figure [TJ] shows the distribution of near-IR luminos- 
ity (from the IRAC 3.6 /im data) with UV luminosity 
for 630 and 344 spectroscopically confirmed galaxies be- 
tween redshifts 1.4 < z < 2.7 and 2.7 < z < 3.7, respec- 
tively, with 3cr upper limits indicated for those galaxies 
that are undetected at 3.6 /im. The effect of the decreas- 
ing fraction of IRAC-detected galaxies with fainter UV 
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Figure 14. Rest-frame 1.1 and 0.9 /im absolute magnitudes (from 
IRAC channel 1 data) vs. UV absolute magnitude for 630 spectro- 
scopically confirmed galaxies at redshifts 1.4 < z < 2.7 (left) and 
344 spectroscopically confirmed galaxies at redshifts 2.7 < z < 3.7 
(right). Points and downward-pointing arrows denote galaxies de- 
tected and undetected (to 3<r), respectively, at 3.6 /im, and the solid 
lines denote the stacked magnitudes as a function of UV luminos- 
ity (Table [8}. The dashed lines show the best-fit linear relations 
obtained from the Buckley-James estimator, taking into account 
both IRAC detections and nondetections. 

luminosities (Figure I15p is evident when examining the 
stacked IRAC fluxes, which are correspondingly fainter 
for UV faint galaxies. The procedure used to compute 
the stacked IRAC fluxes is described in Appendix [DJ 
The trend between stacked IRAC flux and UV magni- 
tude is essentially identical to the trend inferred using 
the Buckley- James estimator on the individual detec- 
tions and non-detections. The best-fit trends between 
near-IR and UV magnitude, taking into account IRAC 
non-detections are: 

Mi.xO ~ 2.3) = (0.90 ± 0.08)M UV - (3.41 ± 1.57), (5) 

with an rms scatter of 1.01 dex, and 



M Q . 9 (z 



> 3.0) = (0.86 ± 0.12)M UV - (3.62 ± 2.45), 



(6) 



with an rms scatter of 1.08 dex. 

As discussed in Appendix [C] the conversion between 
near-IR magnitude and stellar mass depends on UV lu- 
minosity, simply because current star formation can con- 
tribute to the near-IR magnitude. To properly account 
for this effect, we converted the near-IR magnitude found 
for each bin of UV luminosity to a stellar mass using the 
relation between near-IR magnitude and stellar mass ap- 
propriate for that bin of UV luminosity (see Appendix [Ul 
Figure [28l and Table [7]). In this way, we are able to ac- 
count for the star formation history of the average galaxy 
at a given UV luminosity in estimating its stellar mass. 

Figure [16] shows the resulting median stellar masses of 
galaxies in different UV magnitude bins. For compari- 
son, the dashed lines show the results at z ~ 2.30 and 
z ~ 3.05 if we assume a single best-fit linear relation 
between M* and Mi.i,o.9 for all galaxies (irrespective of 
UV luminosity). Interestingly, it is primarily for galax- 
ies brighter than L^ v that the SFR is significant enough 
to bias the M* versus near-IR relation towards higher 
masses. Below L*, the discrepancy in stellar masses de- 
rived using relation for all galaxies versus that derived 
using the relation only for UV faint galaxies, is small. 
This can be attributed to the fact that the former rela- 
tion is dominated by UV faint galaxies at fainter near- 
IR magnitudes (Appendix [C] and Figure [57]) . The rela- 
tionship betwee n UV luminosi ty and mass is similar to 
that derived in iSawickil (|2011| ) for a sample of BX se- 
lected galaxies in the Hubble Deep Field (HDF) North, 
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Figure 15. Fraction of galaxies detected in IRAC channel 1 
(3.6 /xm) as a function of absolute magnitude at 1700 A, at red- 
shifts 1.4 < z < 2.7 (blue) and 2.7 < z < 3.7 (red). Errorbars 
reflect Poisson statistics. 
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Figure 16. Median stellar mass, inferred from stacking the IRAC 
data, in bins of UV absolute magnitude (points and solid lines). 
The stacked IRAC magnitudes were converted to stellar mass as- 
suming the best-fit linear relations between t hes e two quantities in 
bins of UV magnitude (Appendix ICl Figure [28] Table 0. Results 
are also shown (dashed lines) where we have assumed a single best- 
fit linear relation for all galaxies in our sample. The vertical lines 
denote the val ues of L* of the UV lum inosity functions at z ~ 2.30 
and z ~ 3.05 (Rcddy & Stcidcl 2009). The open stars denote the 
median stellar mass in the bins of UV faint photometrically se- 
lected galaxies. The conversion between near-IR magnitude and 
stellar mass for these faint bins is discussed in Appendix [C] The 
dotted black lines indicate lines of constant M/L ratio. 

once we have taken into account differences in SED fit- 
ting by remodeling our ga laxies using the same templates 
(|Bruzual fc Char lotl [20031 ) used in that study. 

The results summarized in Figures [H] and \W\ imply 
that even within the limited dynamic range of UV lu- 
minosity probed with the spectroscopic sample, there is 
a steep trend of mean stellar mass with UV luminos- 
ity. The UV faintest galaxies in our spectroscopic sam- 
ple (around 0.16L*) have mean stellar masses that are 
at least a factor of w 40 smaller than those measured 
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Table 5 

Stacked IRAC Magnitudes for UV faint Sample 



Redshift Interval M(1700A) Range a 3.6 /W 5 4.5 (im b 

1.4<z<2.7 -19.03 -18.53 -19.97 ± 0.10 (393) -20.09 ± 0.09 (260) 

-18.53 -18.03 -19.28 ± 0.26 (170) -19.24 ± 0.34 (112) 

2.7<2< 3.7 -20.05 -19.55 -20.54 ± 0.08 (386) -20.64 ± 0.06 (251) 

-19.55 -19.05 -20.29 ± 0.13 (230) -20.08 ± 0.32 (132) 

a Absolute magnitude range assuming a mean redshift of z = 2.30 and z = 3.05. 
b Uncertainties in absolute magnitude reflect the stacked flux measurement uncertainty 
combined in quadrature with the dispersion in absolute magnitude given the range of 
redshifts of objects in each bin. Parentheses indicate the number of galaxies in the 
stack. 



for UV-bright galaxies with M1700 — —22.5. This steep 
trend is only evident once galaxies undetected in the 
near-IR are included (Figure IT4|): the trend is much shal- 
lower or non-existent when examinin g only those galaxies 
detec ted in if s -band or with IRAC (jShaplev et al.ll200lL 
2005). Similar trends between UV luminosity and stel- 
lar mass have be en found for dropout selected galax ies at 
higher redshifts (jStark et al.ll2009t iLee et al.ll2011al) . and 
could be inferred from the fact that the SFR correlates 
with both UV luminosity (Appendix |B|) and stellar mass 
(Figure [HI). 

7.2.2. Stellar Masses of Photometrically Selected UV Faint 

Galaxies 

To extend these spectroscopic results to fainter UV 
luminosities, we have stacked the IRAC data for pho- 
tometrically selected BXs and LBGs with 1Z > 25.5, 
specifically in two bins with 25.5 < 1Z < 26.5 and 
26.5 < 1Z < 27.0, using the same stacking method dis- 
cussed in Appendix |Dj These photometrically selected 
faint galaxies are likely to lie at the redshifts of in- 
terest, given the strong correlation between contami- 
nation fraction and UV mag nitude (jReddv et all 120081 : 
iReddv fc Steide]|l2009f ). Table[5]lists the stacked 3.6 and 
4.5 /im fluxes for these faint bins, along with the num- 
ber of objects contributing to the stacks. The relatively 
faint stacked IRAC magnitudes obtained for these galax- 
ies suggests that low redshift contaminants, if they exist, 
do not dominate the signal; otherwise, we would have 
expected the stacked IRAC flux to be brighter than the 
observed value. The near-IR magnitudes are converted 
to stellar masses as described in Appendix [Cj These faint 
stacks imply that the trend between stellar mass and 
UV luminosity derived from the spectroscopic sample ex- 
tends to fainter UV magnitudes and lower masses (Fig- 
ure HBJ. 

Having computed their mean stellar masses, we ex- 
amine the UV faint galaxies in the context of the SFR- 
M* relation. The unattenuated SFRs implied by the 
mean UV luminosities of the two bins of photomet- 
rically selected faint galaxies at z ~ 2.3 are ps 2 
and lMoyr" 1 . The average UV slopes for these two 
bins, base d on the relationship between /3 and Mtjv at 
z ~ 2.5 (|Bouwens et all 1201 lal) . are (p) » -1.8 and 
— 1.9. Conver t ing th ese (3 to dust obscurations with the 
iMeurer et al.l (|1999l ) relation and applying them to the 
unobscured (UV) SFRs then implies bolometric SFRs 
of w 5 and 3MQyr _1 . The dispersion in f3 for a given 
Mtjv results in a factor of ~ 2 uncer tainty in these SFRs . 
Adding the intrinsic scatter in the IMeurer et al.l (|1999l ) 



relation increases the SFR uncertainties to a factor of 
w 2.3. Figure IT21 includes the values of inferred bolomet- 
ric SFR for the faint samples analyzed here. 

Given the biases discussed in Sections 16. 4[ O and Ap- 
pendices |B] and [C] it is prudent to determine whether the 
location of the UV faint galaxies on the SFR-Af* plane 
may be biased with respect to all galaxies that lie within 
the same bins of stellar mass. We investigate such biases 
using simulations that are presented in Appendix |El The 
first main result of our simulations is that Malmquist bias 
can result in overestimated mean SFRs in bins of stellar 
mass. One may circumvent this bias by stacking in bins 
of SFR (or UV luminosity) , at the expense of probing the 
intrinsic SFR-Af* relation over a narrower range of M* 
where the flux-limited sample is complete. The critical 
point is that the median stellar mass in bins of SFR may 
not directly translate to the median SFR in that same 
bin of stellar mass. Secondly, in an ideal situation, it is 
desirable to perform full SED fitting for all galaxies in a 
sample, in order to better constrain stellar masses that 
take into account the star formation history and current 
SFR. In our case, a substantial fraction of galaxies in our 
sample are faint (and low mass) and are undetected long- 
ward of the Balmer and 4000 A breaks. One can estimate 
stellar masses for such objects by assuming some M/L 
ratio. However, as we have shown, a proper treatment 
must take into account the UV luminosity dependence of 
the near-IR M/L ratio when estimating stellar masses for 
undetected objects (or objects detected in stack). Hav- 
ing accounted for these biases, we find that, within the 
uncertainties of the faint stacks, the SFR-M* relation ex- 
hibits a close to unity slope to stellar masses as low as 
M* ~ 5 x 10 8 M Q . Progress in quantifying any possible 
evolution of the SFR-Af» relation at low masses should 
be made with future spectroscopic observations of larger 
samples of UV faint galaxies, combined with information 
from deep near-IR selected samples. 

8. DISCUSSION 

A principal aspect of our analysis is the comparison 
of SED-inferred SFRs to those determined from direct 
measurements of dusty star formation at z ~ 2. An im- 
portant conclusion of our analysis is that exponentially 
declining star formation histories yield SFRs that are 
inconsistent with those obtained by combining mid-IR 
and UV data. While both the constant and rising star 
formation histories yield SFRs that consistent with the 
independent measurements, as we discuss below, there is 
additional evidence that suggests that on average galax- 
ies at z > 2 have SFRs that may be increasing with time. 



18 



1000 



100 



> 
D 
+ 
03 



10 



0.1 



I I I i- 

(0(z~2.3)) w =2.4±O.l Gyr" 1 



=L — o t • to i •: • 



10 8 



10 9 



10 10 10" 
M, (M e ) 



10 12 



Figure 17. Specific SFR, </> = SFR/iW*, as a function of stellar 
mass for the MIPS sample of galaxies with redshifts 1.5 < z < 2.6. 
Detections and non-detections are indicated by the filled symbols 
and downward arrows, respectively. The mean specific SFR in bins 
of stellar mass for the spectroscopic sample and photometrically 
selected UV faint sample are indicated by the small open pentagons 
and large open stars, respectively. The weighted mean specific SFR 
over all stellar masses is (<j>) = 2.4±0.1 Gyr" 1 , denoted by the solid 
line. 

8.1. The Typical Star Formation History of 
High-redshift Galaxies 

The SED determined SFRs, and those computed from 
combining mid-IR and UV data, vary roughly linearly 
with stellar mass (Figure H2|) . Hence, as shown in Fig- 
ure[l7l the specific star formation </> = SFR/Af» of galax- 
ies with 1.5 < z S pec < 2.6 is roughly constant over » 2.4 
orders of magnitude in stellar mass up to log[M*/M Q ] = 
11.0, with a weighted mean of <j)(z ~ 2) = 2.4±0.1 Gyr -1 , 
as measured from SFRs[IR+UV]. SFRs[SED] for galax- 
ies with 2.7 < z < 3.7 imply <f>(z ~ 3) = 2.3 ± 0.1 Gyr" 1 . 
Similar values of 4> are obtained if we assume a CSF 
model for deriving the SFRs[SED] and/or stellar masses. 

Several determinations of the specific SFR at low and 
high redshift, nominally derived at a fixed stellar mass 
of 5 x 10 9 M©, are compiled in Figure [HE3 The average 
specific SFRs derived at z ~ 2 and z ~ 3 are similar 
to those measured at redshifts z > 3 (e.g.. | Daddi et al.1 
120071 iStark et all 120091 : IGonzalez et alJl2010t l: including 
dust corrections to these higher redshift points results in 
a roug hly constant specific SFR at z > 4 ()Bouwens et al.l 
l2011al) . and one which is about a factor of 2 — 3 larger 
than the values we find at z ~ 2 — 3. A note of caution 
regarding these higher redshift (z > 3) results is that 
the relatively young age of the universe implies a limited 
dynamic range in M/L ratio, translating to a limited 
range in the possible simple star formation histories. As 
such, it is perhaps not at all surprising to find SFR cx 
M*, given that the young stars that dominate the UV 
continuum also contribute significantly to the rest-optical 
flux, the latter of which is used to constrain the stellar 
mass (i.e., the "outshining" problem as discussed earlier). 

16 The difference in the mean stellar mass, a nd hence mean spe- 
cific S FR, that results from using CB11 versus [Bruzual & CharloU 
(120031') models (the l atter have been used in other studies; e.g., 
IGonzalez et al.1 20101 is negligible compared to the intrinsic dis- 
persion in specific SFR. 
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Figure 18. Evolution of specific SFR <j> as function of redshift for 
galaxies of a stellar mass M, ~ 5 X 10 9 M p). Data are taken from 
the compilation of Gonzalez et al. (2010), which includes values 
fro mlNoeske et ail 1200 7). Daddi ct al. ( 200TD . IStark et ail 120091) ■ 
and Gonzalez et al. (20 10h. We also show the dust corrected values 
from lBouwens et al. (2011a). The rms scatte r in the specific SFR 
at a given mass is take n to be 0.3 d e x for the Epcskc ct al. (120071 ) 
points, 0.32 dex for the lDaddi et ail (120071 ) point, ss 0.3 dex for th e 
dust-corrected values at 4.0 < z < 6 fr om Bouwcns ct al. (2011a), 
and Ri 0.10 dex for the IGonzalez et ail H2010D point at z ~ 7. For 
our determination at z ~ 2.3, we compute a scatter of 0.37 dex 
from a survival analysis that includes MIPS 24 fim non-detections 
(see text; note that this is the scatter of individual galaxies about 
the mean specific SFR, not the error in the mean specific SFR) . A 
similar value of the scatter is adopted for the z ~ 3 determination. 

This degeneracy may also lead to an artificial tightening 
of the scatter in the SFR-Af* relation at high redshift. 
The advantage of the method employed here, for z ~ 2 
galaxies, is that we have (a) determined SFRs largely 
independent of the SED modeling (and hence M#) by 
incorporating the MIPS 24 [im data, and (b) the stellar 
masses are constrained with IRAC data that probe the 
peak of the stellar emission at 1.6 /uru. 

While independent estimates of SFRs are not practical 
at z > 3 given the general faintness of galaxies and their 
expectation of having lower dust content, we note that: 

(a) the mean specific SFR at z ~ 2 is similar to that 
derived at higher redshift, despite the different meth- 
ods used to measure SFRs at z ~ 2 versus z > 3 and 

(b) there is w ell-defined correlatio n between SFR and 
M* as z < 2 (jNoeske et al.ll2007D . where SED model- 
ing independent measurements of SFRs can be obtained 
and where rest-frame near-IR constraints on stellar mass 
are accessible. By extension, these lower redshift results 
hint that there may indeed be an intrinsic correlation be- 
tween SFR and M* at the highest redshifts. Regardless, 
we note that if the stellar masses are severely biased for 
large numbers of galaxies at z > 3 (due to "contami- 
nation" of the rest-optical light from the most massive 
stars), then correcting for this effect result in stellar mass 
densities at z > 4 that would be inconsistent with the 
integrated star for mation history inferred at these early 
cosmic times fe.g.. lReddv fc Stcidel 2009; Bouwcn s et al.1 
I2011al) . Such an inconsistency would be made even worse 
if we then also corrected the stellar mass determinations 
at z > 4 for the effect of strong emi ssion line contribution 
to the bro adband optical fluxes (|Schaerer fc de Barrosl 
120091 120101 ) . With this in mind, we find that even with 
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the iBouwens et all (|2011aD dust corrections, the overall 
picture does not change significantly and the specific star 
formation rates appear to have remained approximately 
constant, within a factor of a few, over the 2.5 Gyr inter- 
val from z — 7 to 2. As such, the average SFRs would be 
roughly proportional to stellar mass, implying that both 
increase exponentially with time. 

8.1.1. Expectations from a Constant Star Formation History 

Given this suggestion, it is useful to place the z ~ 2 — 3 
galaxies in the context of higher redshift samples by 
evolving their SFRs backward in time. For a CSF his- 
tory, the specific SFR will be larger at earlier times. 
For example, a typical galaxy at redshift z = 2.3 with 
cj) « 2.4 Gyr" 1 would have had <j> > 30 Gyr" 1 at z > 2.65, 
w 390 Myr earlier. The upper redshift z = 2.65 is close 
the upper boundary where the 24 ^m data are still sen- 
sitive to the rest-frame 8 /im emission, and the difference 
in time between z = 2.3 and z — 2.65 of « 390 Myr 
is similar to the mean "age" determined for our spec- 
troscopic sample of z ~ 2 galaxies (Figure [5]). Extend- 
ing to higher redshift, the specific star formation would 
increase by more than an order of magnitude in the 
» 80 Myr between z = 2.60 and z — 2.68, reaching 
a value of <f> » 180 Gyr -1 at this latter redshift (Fig- 
ure [T9j) ■ Our sample, which encompasses this higher red- 
shift (z = 2.68), includes very few galaxies with <f> as 
large as the one predicted with a CSF history. The lack 
of galaxies with such high specific SFRs is unlikely to be a 
selection effect because the SFR would remain unchanged 
and such galaxies should still be selectable via their UV 
emission. If anything, the higher redshift progenitors 
of galaxies at z ~ 2 would be more easily detectable 
via their UV emission at the same SFR given that the 
extinction per un it bolometric SFR decreases with in- 
creasing redshift (jReddv et al.l l2010h . Hence, there ap- 
pears to be a true paucity of galaxies with specific SFRs 
4> > 100 Gyr -1 relative to what we would have predicted 
if the majority of galaxies at z = 2.3 had CSF (or declin- 
ing]) histories. Stated another way, the relatively small 
scatter in the SFR at a given stellar mass over the entire 
redshift range 1.5 < z < 2.6 precludes the possibility that 
there exists a large population of galaxies with specific 
SFRs that are substantially larger than the mean values 
observed at these epochs These mean values of <f> are 

17 Aside from the disagreement between SED[IR+UV] and 
SED[SFR] with a declining star formation history, such histo- 
ries also require larger SFRs at earlier times. The tight corre- 
lation between SFR and dust attenuation would then imply an 
anti-correlation between dust attenuation and stellar mass for any 
individual galaxy. If such an anti-correlation existed, we would 
have expected to see a significant evolution in the SFR-Mj rela- 
tion at z > 2, yet thi s is not observed (e.g.. lG~onzalez et al.ll20T(il ; 
IBouwens et al.lf20TTaD . 

18 It is still possible that t here may exist very low m ass galaxies 
with high specific SFRs (e.g.. Ivan der Wei et al. 2011) that do not 
enter our selection, but may eventually merge to produce those 
galaxies seen in our sample. On the other hand, there may also 
exist a heretofore undetected population of low mass galaxies with 
low specific SFRs that merge in the same way as the high specific 
SFR galaxies. In this case, we might expect that the mean specific 
SFR in a given bin of low stellar mass may not be substantially 
different than that of a bin of higher stellar mass. Our IRAC 
stacking analysis supports this expectation to a mass of M, ~ 
5x 10 s Mq. Undoubtedly, future observations should help to clarify 
the situation for galaxies with masses lower than our current mass 
threshold. 
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Figure 19. Same as Figure fl7l where median stacks of the spe- 
cific SFR are shown by the grey points (individual points have been 
suppressed for clarity). The connected red points show the evolu- 
tion in <j> and Af» of a galaxy with A/* = 10 9 , 10 10 , and 10 11 M Q 
and <f> = 2.4 Gyr -1 at z = 2.3, assuming a CSF history. 

roughly independent of stellar mass with (</>} « 2.4 Gyr -1 
(i.e., the lower mass progenitors of z — 2.3 galaxies retain 
roughly the same specific SFR, contrary to our expecta- 
tion if galaxies were evolving with a constant SFR). 

8.1.2. Expectations from a Rising Star Formation History 
and Implications for the Evolution of -Mtjv with 
Redshift 

As discussed in Section 16.21 a rising star formation 
history where the specific SFR <\> is roughly constant 
implies "ages" that can be arbitrarily old if we vary 
the normalization of the star formation history. Given 
such a star formation history, we can estimate the red- 
shift evolution in UV luminosity for each of the galax- 
ies in our sample. For clarity, we show this evolution 
for a subset of 45 galaxies drawn randomly from our 
SED sample of 1959 galaxies, assuming ages that are 
either 2 Gyr or 2.5 Gyr, depending on the redshift of the 
galaxy (so as not to violate the age of the universe at 
each redshift), in Figure [30] The observed UV mag- 
nitudes are calculated assuming that the dust extinc- 
tion varies with total SFR (or intrinsic UV luminosity; 
e.g. iReddv et al.ll2008l;[Bouwens et aD200alReddy et al.1 
120101 ). We adopted the relation between UV luminosity 
and bolometric luminosity in order to apply this correc- 
ti on (see Appendi x jBl adopting the IBouwens et afll2009l 
or IBouwens et al.ll2011al relations between UV slope, fj, 
and Mtjv at z ~ 2.5 yields similar results). The cor- 
rection is renormalized to yield the same dust reddening 
as indicated by the best-fit value of E(B — V) from the 
SED fittingly For comparison, we also show the maxi- 
mum possible UV luminosity (i.e., assuming no dust red- 
dening) as a function of redshift in the left panel of Fig- 
ure 1201 Hence, the curves shown in the left and right 

19 There is evidence of a mild evolution in /3 at a given Mtjv as a 
function of redshift, in the sense that higher redshift galaxies have 
slightly bluer UV slopes than lower red shift galaxies with the same 
absolute UV magnitude (Bouwcns et al. 20lla). For simplicity, we 
have not corrected for this redshift dependence. Doing so would 
shift the UV magnitudes brighter at ear lier times relative to what 
is shown in the right panel of Figure l20l 
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Figure 20. Intrinsic (left) and observed (right) absolute UV magnitude as a function of redshift for 45 galaxies chosen randomly from 
our SED sample. The curves assume a star formation history where the specific SFR <j> is constant. The intrinsic UV magnitude is derived 
directly from the bolometric SFR without any correction for dust. The observed UV magnitude is comput ed under the assumption that 
the dust correction depends on the observed UV magnitude (Rcddy & Stcidcl 2009; Bouwcns ct al. 2011a). The rectangles and red line 
indicate the values and uncertainties in the characteristic absolute magnitude Mfj v derived from the z ~ 3.8 — 8.0 luminosity functions of 
Bouwcns ct al. (2011b) assuming Sz rj 1; and the z ~ 2 — 3 luminosity functions of Rcddv & Stcidcl (2009) assuming Sz = 0.8 and 0.7 at 
2 ~ 2.30 and z ~ 3.05, respectively. 

panels of Figure should encompass the full range of 
possible Mtjy versus redshift tracks for each object. Also 
shown are the characteristic magnitude s M* nr of the UV 
lumin osity functions at z ~ 2 — 3 from IReddv fc Steidell 
(|2009ft and z ~ 4 - 8 from lBouwens et all (|2011b[) . 

In the case of rising SFRs, the redshift evolution in 
Mtjv for z ~ 2 — 3 galaxies with 1Z < 25.5 implies that a 
substantial fraction of them would have been fainter than 
L* (i.e., having sub-L* luminosities) at z > 40 Hence, 
the progenitors of z ~ 2 — 3 galaxies are likely to have ap- 
peared as sub-L* galaxies in dropout samples at z > 4. 
Further, the mean specific SFR of ((f)) ~ 2.4 Gyr -1 at 



z ~ 2 implies a mass doubling time of ~ 300 Myr (Equa- 
tion HJ. Hence, a typical star-forming galaxy at z = 2.3 
has doubled its stellar mass since z = 2.8. Similarly, a 
z = 2.3 galaxy would have had a fifth of its current stellar 
mass (and an intrinsic UV luminosity that is w 1.25 mag 
fainter) at z = 3.0 (e.g., Figures [16] and [20]) . 

At face value, the average rising star formation history 
characteristic of typical star-forming galaxies at z ~ 2 
would imply a significant evolution in L* of the UV lu- 
minosity function at higher redshifts. While this is true 
above z ~ 4, in the sense that L* becomes progressively 
brigh ter with decreasing re dshift from z ~ 7 to z ~ 4 
(e.g., iBouwens et al.lr2011bf ). recent studies suggest lit- 
tle evolution in L* ov er the redshift range z ~ 2 — 4 
IjR.eddv fc Steide]ll2009h . How can we reconcile the rising 
star formation histories of z ~ 2 galaxies with the fact 
that L* does not evolve significantly between z ~ 2 — 4? 
There are two possible explanations. First, observa- 
tions indicate the emergence of a population of mas- 
sive and passively evolvin g galaxies over these redshifts 
(van Dokk um et al.l [20061 ) . implying that some fraction 
of UV-bright gala xies eventually "turn off" by z ~ 2 
(IStark et al.ll2009l) . 

Second, there is a non-monotonic relationship between 
bolometric and UV luminosity for z ~ 2 galaxies, in the 
sense that as the bolometric SFR of galaxy increases, 
dust attenuation will also increase, and eventually there 

20 This conclusion is even stronger for those galaxies that may 
have undergone merging, since the individual components would 
each have an SFR, and hence UV luminosity, lower than that of 
the merged galaxy. 



is a point at which the observed UV luminosity saturates 
IReddv et al.l ((2010). In effect, dust obscuration limits the 
maximum observable UV luminosity of galaxy which, for 
z ~ 2 galaxies, occurs around the value of L* . There- 
fore, as a galaxy's SFR increases with time, there comes 
a point at which dust attenuation becomes large enough 
that its observed UV luminosity no longer increases with 
increasing SFR. It is likely that some combination of 
aforementioned effects (i.e., the fading of some fraction 
of UV-bright galaxies between z ~ 4 and z ~ 2, and the 
limit to the observ ed UV luminosity du e to dust attenu- 
ation; e.g. see also IBouwens et aT]|2009() conspire to pro- 
duce a relatively constant value of L* between z ~ 2 — 4, 
even though the average star formation history may be 
rising during these epochs. 

8.2. Implications for Cold Gas Accretion at High 
Redshift 

In IReddv et al.l (|2006bD . we examined the relationship 
between specific SFR, </>, and cold gas mass fraction, /i: 



,1.4 



4> = C 



i-m' 



(7) 



where we assu med the Kennicutt-Schmidt relation 
(|Kennicuttlll998l) between SFR surface density and cold 
gas mass surface density. The multiplicative factor C in 
the equation above depends on the constant of propor- 
tionality between SFR and cold gas mass surface density 
and the size of the star-forming region. We have defined 
the cold gas mass fraction /i as 



(8) 



We revisit this relationship in light of our new deter- 
min ations of th e spec ific SFRs of z ~ 2 — 3 galaxies. 
The iKennicuttl (|1998l ) relation between SFR and cold 
gas mass surface density reduces to 



M, ~ 5.2 x 10 s 





N 


>l/2 _ 


[M yr-i_ 




kpc 



(9) 



where M g is the cold gas mass is solar masses, ^> is 
the SFR, r±/2 is the half-light radius characterizing the 
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Figure 21. Left: Redshift evolution in gas (solid lines) and stellar (long dashed lines) masses for galaxies with a final stellar mass of 
log M* = 9.5 and 10.5 at z ~ 2.3, assuming a constant specific SFR of <j> = 2.4 Gyr -1 . Also indicated (short dashed lines) is the evolution 



in gas mass fraction for these two examples. Right: Same as left panel where we have assumed 



2.4 Gyr 1 at z < 3.0 based on our 



sample, and ~ 5.0 Gyr 1 at z > 4 based on the dust-corrected specific SFRs presented by Bouwcns ct al. (2011a). We have used a spline 
fit to interpolate the specific SFRs over the redshift range 3.0 < z < 4.0. 



size of the star-forming region, and the power law index 
of the Kennicutt-Schmidt relation is N = 1.4 (below, 
we exp lore how changing the value of N affects our re- 
sults). lLaw et all (|2012D have quantified the sizes and 
their redshift evolution of a subset of the UV selected 
galaxies in the SED sample that have £?5T/WFC3 IR 
imaging. These authors find characteristic half-light radii 
of ri/2 ~ 1.3 and 1.8 kpc for z ~ 2 galaxies in stellar 
mass bins centered at log(M*/Mo) = 9.5 and 10.5, re- 
spectively. These sizes are measured in the rest-frame 
optical and are roughly 20% smaller than those measured 
in the rest-UV, where the latter may be mor e appropri- 
ate fo r the size of the star-forming regions. lLaw et all 
(2012) also find ri/2 oc (1 + z) 1 where 7 ps —1.07, mea- 
sured from z = 3.6 to z = 1.5. This size scaling with a 
power law close to unity is similar to that i nferred from 
highe r redshift (z > 4) dropout samples (Oesc rTet alj 
2010), and is indistinguishable within the errors from 
the rx/2 oc H(z)~ 1 scaling proposed by iFerguson et all 
(2004). In the subsequen t discussion, we assume the re- 
sults of lLaw eF al. (2012) for the redshift evolution in the 
characteristic half-light radii measured in the rest-frame 
UV. 

Adopting Equation [9] we show in Figure [21] the evo- 
lution of the cold gas masses for a typical star-forming 
galaxy with a constant specific SFR of cf> = 2.4 Gyr 
(and, hence, an exponentially rising star formation his- 
tory with T r = <$T X ~ 400MyiE3) and (final) stellar mass 
at z ~ 2.3 of log(M*= 2 - 3 /M Q ) = 9.5 and 10.5 (denoted 
by the solid lines in the left panel of Figure f2l)) . Also 
shown (in the right panel of Figure l2"Tj) is the inferred evo- 
lution assuming a median specific SFR that evolves with 
redshift: at z ~ 2 — 3 we assume 4> ~ 2.4Gyr~ 1 , based on 
our sample; at z > 4 we assume <j6 ~ 5.0 Gyr" 1 , based on 
the dus t-corrected specific SFRs found bv lBouwens et al.l 
(|2011a| ). We have adopted a spline fit to these val- 



21 This exponentially rising timescale is similar to the r r ~ 
0.5 Gyr inferred at z > 5 using the m odels of Bouche etaLI (I2010D ; 
and is also similar to that inferred by Papovich ct al. (2011) to de- 
scribe the average star formation history of galaxies with a constant 
comoving number density of 2 X 10~ 4 Mpc~ 3 at z > 4. 



ues in order to infer the specific SFR over the range 
3-0 < z < 4.0. In addition, the stellar masses are evolved 
backwards in time assuming a gas recycling rate from 
stellar remnants as computed from the CB11 stellar pop- 
ulation synthesis models for an exponentially rising star 
formation history with r r ~ 400 Myr. 

There are several important implications of these re- 
sults. Qualitatively, the inferred evolution assuming a 
constant and evolving median specific SFR (left and right 
panels, respectively, of Figure 12 1 1) are similar. If we 
assume that the Kennicutt-Schmidt relation applies at 
high redshift, then a rising SFR must be accompanied 
by an increase in cold gas surface density. Hence, the 

cold gas mass should increase as » M^ 7 , assuming a 
Schmidt law index of N = 1.4. The cold gas mass frac- 
tion, n, will therefore evolve with time more slowly than 
if the SFR is constant or declining exponentially. The 
results of Figure [21] indicate that the characteristic cold 
gas mass fraction (assuming cf> = 2.4 Gyr -1 ) is relatively 
constant above z ~ 4, and declines more rapidly at lower 
redshift. The average stellar mass of galaxies in our 
spectroscopically confirmed sample at 1.4 < z < 2.7 is 
(M*) = 1.6 x 1O 1O M (see Section [6T]). For this stel- 
lar mass, the cold gas fraction is inferred to be between 
50% - 60% at z > 4 assuming 4> = 2.4 G yr" 1 (or > 60% 
at z > 4 assuming the dust corrected <f> of lBouwens et al] 
I2011al) . and decline to w 40% - 50% atzfs 2.8. More 
generally, on an object by object basis, we find a large 
range in cold gas fraction at z ~ 2, with typical values 
of /i ~ 0.4 — 0.6. This is quite similar to the average cold 
gas fractions of fj, ~ 0.3 — 0.6 inferred for star-forming 
galaxies with similar stellar masses as represented in our 
spectroscopic sample (Figure [9]) at z ~ 2 based either on 
dust-corrected H a (via the Kennicutt-Sch midt relation; 
Erb et alT 2006b') or CO measurements (e.g.jDaddi et alj 



20101 iTacconi et alJ[20Tot iRiechers et al.1 120101 : see also 
Papovich et al.H201lD . 



As we have just discussed, the mean cold gas fraction 
is likely higher at z ~ 4 than at z ~ 2, and the mean frac- 
tion may not evolve strongly above z ~ 4 for typical star- 
forming galaxies at these higher redshifts (e.g., stronger 
evolution may be expected for galaxies with larger stel- 
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Figure 22. Redshift evolution in the SFR (dashed lines) and net 
cold gas accretion rate (solid lines) for galaxies with final stellar 
masses at z ~ 2.3 Af^= 2 - 3 = 10 s , fO 9 , 10 10 , and IO^Mq. The 
curves assume a specific SFR that evolves from <f> ~ 2.4 Gyr - 1 at 
2~2 — 3 to </> ~ 5.0Gyr — 1 at z > 4 and a gas recycling rate 
from stellar remnants computed from the CB11 stellar population 
synthesis models for a exponentially rising star formation history 
with r r ~ 400 Myr. 

lar masses). To illustrate these issues further, we show 
in Figure 1521 the net cold gas accretion rates (as inferred 
from the specific SFRs) and SFRs for galaxies that have 
M z=2.3 = 1Q 8 ) 10 9 j 10 i0 ) and 10 n M by z - 2.3, adopt- 
ing the evolution in </> that was assumed in the left panel 
of Figure [5T] We computed the net cold gas accretion 
rate as M^°' d , which is related simply to: (a) the amount 
of cold gas consumed to form stars in a time 5t (i.e., 
the SFR); (b) the time derivative of the cold gas mass 
(computed from Equation [3]) ; and (c) the rate at which 
the mass formed in stars is released back into the ISM 
as computed from the CB11 stellar population synthesis 
models (R). For a Salpeter IMF, the CB11 models indi- 
cate that the fraction of gas returned to the ISM to the 
stellar mass reaches a value of « 30% at 2Gyr, assum- 
ing an exponentially rising star formation history with 
r r ~ 400 Myr. Hence, we can write the net cold gas ac- 
cretion rate as 

M a c °> d = M g + * - R (10) 

(see also Equation 8 of lPapovich et al.H2011|) . 

Above z ~ 4, the net cold gas accretion rate outpaces 
the SFR by at least a factor of two, with the difference 
between cold gas accretion rate and SFR being largest for 
lower mass galaxies at any given redshift above z ~ 3 — 4. 
It is evident that the SFR becomes comparable to the net 
cold gas accretion rate at a redshift which is dependent 
upon stellar mass (i.e., the redshift where the dashed and 
solid lines come together in Figure |2"2"|) . In particular, 

*f? « Ma°c d & t higher redshifts for galaxies with larger 
stellar masses. 

Above z ~ 2 — 3, the accretion rate outpaces the SFR 
and the average cold gas fraction r emains consistently 
high (e.g., right pa nel of Figure |2"T1 iBouche et al.l [20101 : 
see also Figure 4 of lPapovich et al.ll2011[ ). Of course, the 
actual dispersion in fJ, at a given redshift (even above 
2^4) may still be quite large. Nonetheless, the steep 



UV l uminosity function at z > 4 (e.g., iBouwens et al.l 
2009) implies that less U V luminous and h e nce lower 
stellar mass galaxies (e.g., iStark et "ail 120091 : iLee et all 
l2011alf bl: see also Figure El) should dominate, and hence 
the number- weighted mean cold gas fraction at these red- 
shifts should be quite high with fi > 60% (Figurel2T|). Fi- 
nally, we note that the difference in M^ d and SFR will 
be slightly l arger at z > 4 h ad we assumed a shallower 
slope of the iSchmidu (119591) relation of N ~ 1.1 - 1.2 
(e.g., as suggested bv lGenzel et al.ll201Q[ ). with a further 

result being that the point at which \1/ ps M^r shifts to 
lower redshifts for galaxies of a given stellar mass. 

The physical picture implied by the redshift evolution 
in cold gas accretion is that on average galaxies must 
be continually accreting gas to support their increasing 
SFRs. Indeed, the gas consumption timescales assuming 
a CSF history for z ~ 2 UV selected galaxies are on 
average a factor of two lower than median stellar ages, 
implying that signi ficant gas accretion must be occurring 
at z ~ 2 (|Erbl2008f ). while at the same time, such galaxies 
are inferred from clustering measurem ents to have long 
duty cycles (e.g., iConrov et al.1 120081 ) . A high rate of 
gas accretion is also required to explain the increasing 
SFRs (and SFR density) cont ributed by galaxies of a 
fixed number density at z > 4 (Pa povich et alJfeOllD . 

We conclude with a caution against over-interpreting 
the physical implications of the redshift at which \P « 
M™1 . One might speculate that below this redshift 
which for typical galaxies with M* > 10 9 M Q occurs 
around z ~ 2 (Figure |2"2"T) — the SFR declines primar- 
ily in response to a lower cold gas accretion rate. While 
this may be the case — and would explain naturally the 
overall decline in the cosmic SFR density at the same 
redshifts (z ~ 2) — as we discuss in the next section, the 
cold gas accretion computed above evidently constitutes 
only a small fraction of the total baryons accreted onto 
halos. Therefore, some other mechanism(s) must be re- 
sponsible for regulating the large number of baryons that 
never end up in stars at high redshift. We discuss this 
further in Section 15731 

8.3. A Consideration of Gas Inflows, Outflows, and the 
Expectations from Simulations 

So far, we have not distinguished between gas inflows 
and outflows, and have only referred to the net cold gas 
accretion rate as inferred from the specific SFRs of z > 2 
galaxies. We can estimate the cold gas inflow rate by 
making some assumptions of the gas outflow rate. Taking 
the net cold gas accretion rate as the difference between 
inflow (Mi n ) and outflow rates (M out ), and assuming that 
Mout — Vout ^ , where r ?o„t is commonly referred to as the 
"mass-loading factor, "0 we can rewrite Equation [TOl 

M in = M g + (1 + 7? QUt )* - R. (11) 

22 Hydrodynamical simulations suggest that under the assump- 
tion of momentum-driven winds, the instantaneous mass outflow 
rate ?7inst °c , where v c is the circular velocity, and hence 
r? inst oc M~ 1/3 (e.g., IFinlator & Dav^l200l ; IDave et al.|[20TlT) . In 
this case, ?7 ou t could be thought of as the mass-lo ading factor aver- 
aged over the star formation history of the galaxy (Finlator & Dave 
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Figure 23. SFR (dashed line) and inferred cold gas inflow rates 
(solid lines) for different values of the mass-loading factor rj, for a 
galaxy with a final stellar mass of M, = 1.6 X 10 10 Mq at z = 2.3. 
The curves assume a specific SFR that evolves from cf> — 2.4 Gyr 
at z ~ 2 — 3 to 4> — 5.0 Gyr -1 at z > 4 and a gas recycling rate 
from stellar remnants computed from the CB11 stellar population 
synthesis models for a exponentially rising star formation history 
with T r ~ 400 Myr. The hashed red region denotes the inferred 
baryonic accretion rate (dark matter accretion times the baryonic 
fraction of fi b /fi m = 0.165) for halos with log[M h /M ] ~ 12.0 - 
12.2 at 2 = 2.3. 

Direct measurements of ?7 ou t in high-rcdshift galaxies are 
sparse, given the difficulty of measuring the covering frac- 
tion and terminal velocity of the gas and its radial extent, 
as well as the HI column density (due to saturation of 
the typi cal absorption lines used to make such measure- 
ments). iPettini et al.l ([20021 ) estimated the mass outflow 
rate of the low-ionization and HI gas of the lensed Ly- 
man Break galaxy cB58, finding 77 ou t ~ 1.75. This is 
likely a lower limit because a significant fraction of the 
outflowing gas may be at higher vel ocity and / o r in a more 
highly ionized state. Alternatively, lErbl (|2008D finds that 
77 out ~ 1 is required to reproduce the mass-metallicity re- 
lation at z ~ 2. From the theoretical standpoint, hydro- 
dynamical simulations find mean mass-loading factors of 
at least a factor of a few for galaxies with stellar masses at 
the low end of those found in our sample (M* ~ 10 9 Mq), 
and being less than unity for the most massive galax- 
ies (lOppenheimer fe Dav| 120061 [20081: iFinlator fe Pavel 
120081: lOppenheimer et aiTl201 It iDave et al.l 1201 ID . To 
encompass the range of observed and expected mass out- 
flow rates, we consider the effect of changing 77out on the 
cold gas inflow rate inferred for a typical galaxy in our 
sample with M^ =2 - 3 ~ 1.6 x 1O 1O M0 , as shown in Fig- 
ure pi 

Figure [23] suggests that the total baryonic accretion 
rate must be quite large relative to the SFR if r? ut 
is at least a factor of a few. It is useful to com- 
pare these inferred inflow rates to the total gas accre- 
tion expected for the dark matter halos that host galax- 
ies in our sample. Based on clustering analysis, the 
correlation length of galaxies in our spectroscopic sam- 
ple (with luminosities around L* of the UV luminos- 
ity function; IReddv fe Ste"idel[2009l) is 6.5 ± 0.5 h~ x Mpc 
(comoving), which corresponds to dark matter halos 
with masses log[Mh/M Q ] > 11.8, with an average halo 
mass of (log[Mh/M Q ]) = 12.2 and median halo mass of 



log[Mh/M ] = 12.0 (Trainor & Steidel, in preparation). 
The total baryonic accretion rate, which is a generic and 
well-determined quantity from ACDM cosmological sim- 
ulations, can be approximated as 



M b ~ 6.6 



My, 



10 12 Mr: 



(l + z) 2 - 25 /o.i65 Moyr- 1 (12) 



(jDekel et al.ll2009t ). where /o.i65 is the fraction of baryons 
to total matter (baryons plus dark matter) in units 
of the cosmological value, fib/^m = 0.165 ± 0.009 
(jKomatsu et all I2009D FI Employing Equation [HI we 
show in Figure [53J the total baryonic accretion rate (red 
shaded region) based on hydrodynamical simulations for 
the halos that are inferred to host galaxies in our sample. 

The expected baryonic accretion rates from the simula- 
tions are relatively constant with redshift, owing to the 
fact that the larger accretion rates expected for higher 
mass halos at lower redshift are counteracted by cosmo- 
logical expansion (which tends to decrease the accretion 
rate at a given halo mass). At z ~ 2 — 3, larger values 
of the mass-loading factor (77 > 1) yield a total cold gas 
inflow rate (as inferred from the specific SFR) that is 
similar to the baryonic gas accretion expected for halos 
with log[M h /M ] ~ 12.0 - 12.2. Strikingly, however, at 
z > 3.5 there is a large discrepancy between the baryonic 
gas accretion rate expected for the halos and the cold gas 
accretion rate inferred from the specific SFR with reason- 
able values of rj (cyan and orange lines in Figure [231) . For 
instance, there is a factor of ~ 40 difference between the 
two expectations at z ~ 6. 

While the discrepancies at z > 3.5 could be resolved 
by adopting a mean mass-loading factor that is red- 
shift dependent and substantially larger than rj — 4.0, 
several other possibilities have been suggested in the 
literature. At face value, the hydrodynamical simula- 
tions suggest that the SFR for a galaxy in a given halo 
should be roughly constant with time if the fraction of 
accreting baryons that turns into stars is roughly con- 
stant with redshift. Could the galaxies in our sample 
truly be forming stars at a constant rate? As argued 
in Section 18.1. 1[ this is unlikely given the fact that we 
do not observe galaxies at the higher redshift end of 
our selection function with specific SFRs that are sub- 
stantially larger than the median specific SFR. On the 
other hand, the specific SFRs measured at z > 4 are 
still subject to several systematic uncertainties, includ- 
ing the overestimation of stellar mass due to emission 
line contribution to the broadba nd near-IR luminosity 
(e.g.. iSchaerer fc de Barrosll2010D . A star formation his- 
tory that rises more slowly than an exponential (i.e., as 
inferred from mean specific SFRs that increase towards 
higher redshift) can bring the inferences from the spe- 
cific SFR closer to the expectations from the simulations. 
Alternatively, an IMF that becomes progressively more 
bottom-light towards higher redshift could also lead to 
elevated specific SFRs at early times, though such IMF 
evolution would result in an integrated SFR density that 

23 The fraction, t, of the baryonic accretion rate that is in a 
cold state (T < 10 4 K) is model dependent and highl y uncertain , 
with suggested values of e ~ 0.4 — 0.8 (e.g., [Bouchc et af] 120101 ; 
Fauchcr-Gigucrc et al. 2011). Adopting these typical values would 
lower the total cold baryonic accretion rate by a factor of < 2 and 
does not change any of the subsequent conclusions. 
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disagr ees with estimates of t he stellar mass density at 
z > 4 (IReddv fc Steidell[2009h 

From the theoretical side, iKrumholz fc Dekel (|2011[ ) 
suggest that star formation at z > 2 may be suppressed 
relative to the baryonic accretion rate due to limited 
free-fall time and a lack of metals in halos of mass 
Mh < 10 11 Mq. This suppression results in a departure 
bet ween the gas accr etion rate and SFR as predicted from 
the lKennicuttl (|1998| ) relation. However, we note that the 
rcdshift and mass dependence of the baryonic accretion 
rate implies that the halo of an L* galaxy at z <~ 2.3 with 
Af h 10 12 M , would have had M h ~ 1.6 x 10 11 M at 
z ~ 6, on average. This mass lies just at the threshold 
of where metallicity effects are infer red (from the simu- 
lations) to suppress star formation (Krumh olz fc Dekel 
1201 If ). Yet, even for halos as massive as Mh ss 10 11 M Q 
at z ~ 6, we infer gas inflow rates that are a factor of 
40 — 200 times smaller than those inferred from the sim- 
ulations (Figure ESjEB 

An other possibility suggested by ([Krumholz fc Dekel 
I2011h is that the bulk of the stellar mass growth at high 
rcdshift occurs from the accretion of stars formed in other 
galaxies, in which case the observed specific SFRs would 
be biased lower than the intrinsic specific SFR derived 
by considering only in situ star formation. Nonetheless, 
the actual SFR of galaxies within such halos apparently 
constitutes only a small fraction of the baryons accreted 
onto those halos. It is possible that some combination 
of the aforementioned effects may result in the very ef- 
ficient suppression of star formation in even relatively 
massive halos (Afh ~ 10 1:l Mq) at z > 2, but we note 
that such mechanisms must become increasingly more 
important with increasing redshift. Irrespective of the 
exact mechanism, what is perhaps quite striking from 
Figure [221 is that star formation in general may have 
been an extremely inefficient process even at early cosmic 
times when galaxies were first assembling. 

The relative inefficiency of star formation at early 
times can be envisioned also by examining the stellar- 
to-halo mass ratio (M*/Mh), as shown in Figure [2"4l At 
z = 2.3, the mean stellar mass of galaxies in our sample, 
and their median and average halo masses as inferred 
from clustering, imply M*/Mh ~ 0.10 — 0.16. Evolving 
the stellar masses backwards in time assuming the typ- 
ical star formation history (e.g., as in Figure 121"]). and 
evolving the halo masses toward higher redshift assum- 
ing the baryonic accretion rate formula (Equation I12[) , 
then implies a stellar-to-halo mass ratio that is an order 
of magnitude smaller at z ~ 4.5 than at z = 2.3, and 
a factor of 40 lower at z ~ 6 than at z — 2.3. The re- 
sults imply again that on average a very small fraction 
the baryons accreted onto halos at early times actually 
ends up forming the stellar mass of a galaxy at simi- 
lar epochsE3 These results agree qualitatively with the 



10 12 M halos at z 
ill - 



It is possible some fraction of the ws ±u~ iviq 
2.3 had progenitor halos that are less massive than Mj, ~ 10 11 Mq 
at 2 ~ 6, where metallicity effects would be noticeable, if such 
low mass halos merge in a way that the time-averaged baryonic 
accretion rat e evolution occurs more rapidly than that assumed by 
Equation 1121 

25 We have been careful to emphasize that only a small frac- 
tion of baryons accreted in the early history of a typical galaxy 
contributes to star formation at those epochs. This does not neces- 
sarily imply that such baryons are forever precluded from forming 
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Figure 24. Stellar-to-halo mass ratio as a function of redshift 
for a typical star-forming galaxy with a stellar mass at z = 2.3 of 
M* =2 3 = 1.6 X 10 10 Mq. Clustering results are used to estimate 
the corresponding average and median dark matter halo masses 
(Trainor & Steidel, in preparation). The stellar masses are evol ved 
backwards in time using the same assumptions as in Figure 1211 
The median and average halo masses are evolved to higher redshift 
assu min g the baryonic rate accretion formula specified by Equa- 
tion [12] 

lower M*/Mh expected at higher redshifts, and for lower 
mass halos at a giyen re dshift, from A^-body simulations 
(e.g.. lMoster et al.|[20T0h . 

8.4. Summary and Implications 

In this section, we have discussed the typical star for- 
mation histories of high-redshift galaxies based on (a) 
an analysis of the median specific SFR as a function of 
stellar mass at z ~ 2, combined with (b) the observa- 
tion that the specific SFR does not evolve strongly with 
redshift above z ~ 4, at least for galaxies for which such 
measurements have been obtained. These two pieces of 
evidence lead us to the conclusion that on average the 
SFRs of z > 2 galaxies are likely rising monotonically 
with time. In this case, the "age" of the galaxy is an 
ill-defined quantity and in principle could be arbitrarily 
old depending on the normalization of the star formation 
history. As such, it is possible that a substantial fraction 
of the UV faint sub-i* galaxies at z > 4 evolve into UV- 
bright galaxies (i.e., ones that we would spectroscopically 
confirm) at z ~ 2. Studies of the clustering properties of 
sub-L* galaxies combined with halo abundance matching 
will undoubtedly help elucidate the probable progenitors 
of z ~ 2 galaxies. 

We proceed to examine the evolution in gas masses 
and gas accretion rates implied by rising star formation 
histories assuming that the Kennicutt-Schmidt relation 
holds at high redshift. The results imply that the gas 
masses will increase with increasing SFR, and that the 
gas mass fraction will evolve more slowly than in a sce- 
nario in which galaxies are forming stars at a constant 
or declining rate. We find in general that mean cold gas 
fraction remains quite high at > 60% at z > 4. A further 
consequence of the slow evolution of gas mass fraction at 
z > 4 is that the average gas phase metallicity may not 

stars at some later epoch (e.g., at z < 2), particularly if the gas is 
continually cycled into and out of galaxies. 
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evolve with time as strongly at z > 4 relative to that ob- 
served at z < 2. This is because newly accreted gas will 
continually dilute the overall gas-phase metallicity. At 
face value , the fundamental me tallicity relation (FMR) 
derived in lMannucci et all ([2010) implies an s» 0.5 dex in- 
crease in mean oxygen abundance from z = 7 to 4. How- 
ever this relation has only been calibrated up to z ~ 2.5 
and for larger SFRs at higher redshifts, than the SFRs 
considered here at z > 4. Direct metallicity measure- 
ments for such high-redshift galaxies, while unobtainable 
at the present time, will prove crucial in connecting the 
metallicity evolution with that expected from the aver- 
age increasing star formation rate between z = 7 and 
2. 

At these early times, the net gas accretion rate ex- 
ceeds the SFR by at least a factor of two, where the 
former is determined by relating the sp ecific SFR to the 
cold gas mass via the IKennicuttl ()1998[ ) relation. For the 
progenitors of typical star-forming galaxies in our z ~ 2 
sample, the SFR eventually equals the net gas accretion 
rate around z ~ 2 — 3, at which point the gas fraction 
decreases below 50%. These results are consistent with 
measurements of the gas mass fractions of high-redshift 
galaxies, and the fact that such galaxies must be accret- 
ing substantial amounts of gas in order to sustain their 
SFRs. 

We proceed to compare the gas accretion rates inferred 
from the specific SFR — combined with some assumption 
for the mass- loading factor of outflowing winds — with the 
accretion rates expected for the dark matter halos that 
are inferred from clustering to host galaxies in our sam- 
ple. We find that the two methods of inferring gas accre- 
tion rates agree at z < 3, but they become increasingly 
divergent above z > 4, such that less than one percent 
of the baryons accreted by a log[Mh/M Q ] = 11.2 halo 
at z ~ 6 would end up in stars at that epoch. While 
several possibilities are discussed to explain this discrep- 
ancy, the results imply that star formation must have 
been extremely inefficient at early cosmic times. 

9. CONCLUSIONS 

We use a large spectroscopic sample of L* galaxies at 
redshifts 1.4 < z < 3.7, and a photometrically selected 
sample of fainter galaxies at the same redshifts, to con- 
strain the average star formation history of typical star- 
forming galaxies between redshift z = 7 and z = 2. Our 
analysis takes advantage of rest-frame UV spectroscopy, 
and optical, ground-based near-IR, Hubble Space Tele- 
scope WFC3/IR, and Spitzer Space Telescope IRAC and 
MIPS 24 (Ltm imaging, in order to measure the stellar pop- 
ulations and dust obscured star formation of galaxies in 
our sample. We have taken into account a number of sys- 
tematic effects in computing these quantities, including 
the contribution of line emission to the broadband pho- 
tometry, dynamical time constraints on galaxy ages, and 
the stellar population dependence of the UV attenuation 
curve. 

With these data, we perform detailed comparisons be- 
tween SED-inferred SFRs (SFRs[SED]) and those com- 
puted by combining deep UV and MIPS 24 fixn data 
(SFRs[IR+UV]), for 302 spectroscopically confirmed 
galaxies at redshifts 1.5 < z < 2.6, where the MIPS 
24 fim observations are sensitive to the dust emission fea- 
tures at rest-frame 8/xm; such comparisons allow us to 



rule out some simple star formation histories. The uncer- 
tainties in stellar masses and ages are discussed. We pro- 
ceed to examine the correlation between SFR[IR+UV] 
and stellar mass, as well as SFR[SED] and stellar mass 
for a larger sample of 1959 spectroscopically confirmed 
galaxies, taking into account several systematic effects in- 
cluding Malmquist bias and differences in the conversion 
from near-IR luminosity to stellar mass. These results 
are extended to lower SFRs and stellar masses via a pho- 
tometrically selected sample of 563 UV faint (TZ > 25.5) 
galaxies over the same redshifts. This information is 
used to constrain the average star formation history of 
high-redshift galaxies, which in turn yields information 
on the gas mass evolution and gas mass accretion rates 
for galaxies of different stellar masses. 
The main conclusions are: 

• SED modeling that assumes exponentially declin- 
ing star formation histories results in SFRs [SED] 
that are on average a factor of 5 — 10 x lower than 
SFRs[IR+UV]. On the other hand, modeling typ- 
ical star-forming galaxies (Lboi 10 12 L Q ; ages 
< 100 Myr) wit h constant or r i sing st ar formation 
histories and a iCalzetti et al.1 (|2000f ) attenuation 
curve results in SFRs [SED] that are in good agree- 
ment with SFRs[IR+UV]. This suggests that most 
L* galaxies at z ~ 2 — 3 did not follow a declining 
star formation history prior to the epoch at which 
they are observed. The traditional use of exponen- 
tially declining models to describe the star forma- 
tion histories of high-redshift galaxies is called into 
question by the analysis presented in this paper. 

• Assuming standard modeling assumptions 
(|Calzetti et al.l |2000| attenuation curve and no 
age constraint) results in SFRs [SED] that are 
systematically larger than SFRs[IR+UV] for 
galaxies inferred to have younger stellar pop- 
ulations (< 100 Myr). However, by building a 
coherent picture for the dust attenuation (SMC- 
like) and dynamical timescale of younger galaxies, 
we are able to reconcile the SFRs [SED] and 
SFRs[IR+UV] for such galaxies. 

• Rising star formation histories yield stellar masses 
that are comparable to (and ages that are typically 
older than) those obtained for a constant star for- 
mation history, owing to the slower evolution in 
M/L ratio for rising histories. More generally, the 
"age" of a galaxy is an ill-defined quantity for ris- 
ing histories, as the observed M/L ratio can be 
obtained by simultaneously varying the age of a 
galaxy and the normalization of its star formation 
history. 

• The scatter in the M/L ratios of galaxies at a 
given luminosity decreases at longer wavelengths. 
Nonetheless, the scatter, which is induced partly by 
the current star formation, is still relatively large 
(0.20 dex) at wavelengths where the stellar emission 
peaks (> 1 /an). The luminosity dependence of the 
M/L ratio at 1.1 /im is partly a result of sampling 
a finite range of UV luminosity in a flux limited 
sample like ours, combined with the effects of dust 



obscuration. The systematic variation of the stel- 
lar mass versus near-IR magnitude relation with 
UV luminosity can result in biased measurements 
of the masses of UV- faint galaxies. 

Taking into account IRAC upper limits and the 
aforementioned UV luminosity dependence of the 
conversion between near-IR magnitude and stellar 
mass, we find a strong correlation between UV lu- 
minosity and stellar mass, even within the limited 
dynamic range of our spectroscopic sample. An 
IRAC stacking of photometrically selected UV faint 
galaxies to 1Z ~ 27.0 shows that this trend extends 
to UV faint galaxies. In particular, 0.1L* galax- 
ies at z ~ 2 have stellar masses that are a factor 
of w 140 x smaller than those of the UV brightest 
(~ 51/*) galaxies in our sample. 

Combining the UV-faint data with our spectro- 
scopic sample, and taking into account systematics 
in the calculation of stellar masses, we find a close 
to unity slope of the SFR-M* relation over 2.5 or- 
ders of magnitude in stellar mass and 2.8 orders 
of magnitude in SFR at z ~ 2. This relation im- 
plies that the specific SFR (SFR/M*) is roughly 
constant for stellar masses M* > 5 x 10 s M Q and 
SFRs > 2M Q yr- 1 . 

The constant specific SFR with mass at z ~ 2, com- 
bined with the fact that this median SFR is similar 
to that found at z > 4, implies that the SFRs must 
be increasing with (and roughly proportional to) 
the stellar mass. A constant star formation history 
would imply much larger specific SFRs at slightly 
higher redshifts, and such objects are not observed 
at those higher redshifts that are probed within our 
sample. If galaxies on average have SFRs that in- 
crease with cosmic time, then it suggests that the 
UV-bright galaxies observed at z ~ 2 — 3 had pro- 
genitors that would be detected as sub— L* galaxies 
in z > 4 dropout samples. 

If we assume tha t the Kennicutt-Schmidt relation 
iKennicuttl (|1998f ) holds at high redshift, then an 
average rising star formation history implies a gas 
mass that increases with time and a gas mass frac- 
tion that evolves slowly at z > 4. The specific SFRs 
at z > 2 indicate a net cold gas accretion rate that 
outpaces the SFR at high redshift, until z ~ 2 — 3, 
when the SFR becomes comparable to the net cold 
gas accretion rate for typical star-forming galax- 
ies. We have inferred the cold gas inflow rate from 
the net cold gas accretion rate with some assump- 
tion for the mass-loading factor of galactic outflows. 
The inflow rates derived in this manner diverge in- 
creasingly from the total gas inflow rates expected 
for the dark matter halos that host our galaxies (as 
inferred from clustering analysis; Trainor & Stei- 
del, in prep.). At face value, the observations im- 
ply that less than one percent of the baryons ac- 
creted onto halos of mass Mh ~ 1.6 x 10 11 M at 
z ~ 6 end up contributing to star formation at that 
epoch, thus highlighting the inefficiency of star for- 
mation even at early cosmic times when galaxies 
were growing rapidly. 



There are two important points to make. First, taken 
together, the results of the SFR comparisons (Section [4}, 
the near-unity slope of the SFR-M* relation, and the 
roughly similar specific SFRs at z > 2, point to a pic- 
ture in which the average SFRs of high-redshift galaxies 
increase with cosmic time. This behavior produces natu- 
rally the ob served rise in the SFR density between z = 7 
and z — 4 (jPapovich et al.|[2011f) . Second, an analysis 
of the specific SFRs indicate that these rising histories 
cannot be sustained indefinitely, and once the SFR ex- 
ceeds the gas accretion rate, we expect the gas mass to 
decrease leading to an eventual decrease in SFR. How- 
ever, our analysis underscores the fact that the net cold 
gas accretion inferred from the specific SFRs is only a 
small constituent of the total baryons accreted onto ha- 
los at z > 2, and that in general star formation must 
have been very inefficient at early times. 

We conclude by noting several future avenues of in- 
vestigation to confirm more robustly our findings. We 
have emphasized a few of the biases (e.g., Malmquist 
bias and systematics in the conversion between near-IR 
luminosity and stellar mass) that are important for anal- 
yses of the SFR-M* (and specific SFR versus M*) re- 
lations. The presence of a Malmquist bias, however, 
does not preclude the possibility of there being a real 
evolution in the SFR-M* relation at faint luminosities 
and low stellar masses. Hence, significant progress in 
evaluating possible evolution can be achieved with deep 
near-IR photometry (e.g., with Hubble Space Telescope's 
WFC3/IR Camera) that in turn enables SED-fitting for 
individual UV-faint and low mass galaxies at high red- 
shift. As emphasized above, there is a substantial disper- 
sion in the star formation histories of individual galax- 
ies; future deeper studies will enable us to quantify how 
the dispersion in SFRs may change as a function of M* , 
particularly for UV-faint galaxies at z > 3 which, from 
clustering studie s, may exhibit a more episodic star for- 
mation histor y (iLee et all 1201 lbl). Further, as recently 
highlighted in lSchaerer fc de Barrosl (|2009ll2010ft . the ex- 
pected high equivalent width emission lines (e.g., Ha) of 
z > 4 galaxies can result in factors of 2 — 3x overesti- 
mates of their stellar masses from broadband SED fitting. 
Correcting for the presence of those emission lines may 
potentially result in a stronger evolution of specific SFR 
than has been found in other studies of z > 4 galaxies 
dStark et alJ 12009b iGonzalez etafl 120101: IBouwens etHI 
1201 laf) . and may partially resolve the discrepancies be- 
tween current observations and numerical simulations re- 
garding the evolution of specific SFRs at z > 4 (e.g. , 
[ BoucM et al.ll2010HDutton et al.H2010b iDaye et alj |2011b 
iWeinmann et alj|2011b iKrumholz fc Dekelll2011t see dis- 
cussion in IBouwens et al.1l2011al ). The forthcoming deep 
near-IR imaging surveys, along with deep imaging and 
spectroscopic campaigns with the next generation of 
ground-based (TMT and GMT) and space-based facili- 
ties (JWST), promise substantial progress in understand- 
ing star formation and the buildup of stellar mass at early 
cosmic times. 
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A. CONVERSION FROM UV LUMINOSITY TO STAR FORMATION RATE 

The conversion from UV luminosity to SFR will depend on the star formation history and age of a galaxy. This is 
illustrated in Figure [25j where we show the ratio of SFR to the 1700 A luminosity as a function of age for different 
sta r formation hi s tories . For ages < lOM yr, the factor needed to convert a UV luminosity to a SFR is larger than 
the iMadau et al.1 (|1998f ): IKennicuttl (|1998t ) values because the ratio of O stars to B stars (both of which contribute 
significantly to the rest- frame 1700 A luminosity) is larger than the equilibrium ratio reached after the B star main 
sequence lifetime of ~ 100 Myr. For older ages and exponentially declining star formation histories, the UV luminosity 
does not fall off as quickly with time as the SFR. In these cases, the ratio of O stars to B stars is lower than the 
case of continuous star formation (i.e., B stars will still contribute significantly to the UV luminosity even after the 
production rate of O stars declines). As the decay timescale Td becomes larger and iago/^d < 1> the conversion factor 
approaches that of the CSF history. Similarly, exponentially rising star formation histories have a conversion factor 
similar to, but slightly larger than, those obtained for a CSF history. In the rising histories, the ratio of O and B 
stars will stabilize to a value larger than in the CSF case since the SFR is continually increasing. The conversion 
factor for rising histories with larger r r approaches the factor for a CSF case. For convenience, we have tabulated the 
form of the simple polynomial fits that describe the relationship between the ratio of SFR to 1700 A luminosity as a 
function of age for the different star formation histories, assuming the latest CB11 models with a Salpeter IMF and 
solar metallicity (Table |6|). In this parameterization, the conversion factor between UV luminosity and SFR is CuV) 
where 

SFR £1700 , . 

tt — = Cuv x (13) 

M yr 1 erg s 1 Hz 1 

and Cuv is approximated with a polynomial function of £ = log(£ agc ): 

5 

log[Cuv(£)] = $>[/(£)]' (14) 

i=0 

where /(£) = £ or l/£, and the coefficients a; are given in Table[6l 

In the present context, changes in the conversion between UV luminosity and SFR will not significantly affect the 
bolometric SFRs for two reasons. First, as discussed in Section [U the range of ages and star formation histories 
that are plausible at z ~ 2 rule out conversion factors that vary by more than w 20% from the equilibrium value. 
Secondly, the bolometric SFRs are less sensitive to changes in the UV-SFR conversion factor given that, for most of 
the gala xies in our sample, a sign ificant fraction of the bolometri c luminosity (~ 8 0%) is emergent in the IR and not 
the UV (jReddv et alj|2006bl l201Cft . For simplicity, we assume the IKennicuttl (|1998t ) relation to convert UV luminosity 
to SFR, and we discuss in Section [4] how changing the conversion factor between these quantities affects (or does not 
significantly affect) our results. 

B. THE EFFECT OF MALMQUIST BIAS ON THE MEAN SFR AT A GIVEN STELLAR MASS 

In this section, we explore the effect of Malmquist bias from a flux-limited sample on the mean SFRs in bins of 
stellar mass. To illustrate this bias, we simulated the effect of imposing a spectroscopic limit on a SFR versus stellar 
mass trend with an intrinsic slope of unity, and an rms dispersion of 0.38 dex (based on the dispersion around the 
best-fit correlation shown in Figure IT2"]). Because our sample is selected based on (unattenuated) UV emission, we first 
quantified the relationship betwee n bolometric SFR and bolometric luminosity using the 302 galaxies in the MIPS 
24 /jm sample. iReddv et al.l (|2010H give the relationship between bolometric and UV luminosity for this sample (with 
a mean redshift of (z) ~ 2) as 

log[£uv/L ] = log[L bol /L ] - log[L bol 1 / a 10- b/a + 1], (15) 

where a = 0.69 ± 0.03 and b — 10.91 ± 0.04. This relationship is redsh ift dependent, because the ratio of bolometric 
to UV luminosity at a given SFR decreases with increasing redshift (jReddv et al.ll2006bl [2010). The mean trend 
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Table 6 

Coefficients of Polynomial Fits to Logarithm of SFR/UV Conversion Factor, log[Cuv] 



SFH a 


7" 

(Myr) 


f(0 h 




Ol 


«2 


«3 


0,4 


0-5 


Age Fiange c 
(Myr) 


Constant 


oo 




-27.822 


-0.969 


2.511 


-2.000 


0.704 


-0.089 


> 3 


Declining 


10 




-26.658 


-1.662 


0.568 


-0.689 


0.969 


-0.395 


3 - 255 


20 


£ 


-26.433 


-3.065 


3.808 


-3.592 


1.936 


-0.425 


3 - 508 




50 


£ 


-25.634 


-7.271 


11.088 


-8.624 


3.237 


-0.471 


3 - 1433 




100 


t 


-26.309 


-3.416 


3.665 


-2.245 


0.716 


-0.093 


3 - 905 




200 


i/e 


-28.906 


3.768 


-5.127 


3.615 


-1.176 


0.140 


3 - 1900 




500 


m 


-28.008 


-0.238 


1.404 


-1.218 


0.449 


-0.059 


3 - 4000 




1000 


ve 


-27.889 


-0.725 


2.162 


-1.764 


0.630 


-0.081 


3 - 7000 




2000 


i/e 


-27.847 


-0.888 


2.409 


-1.939 


0.687 


-0.088 


> 3 




5000 


V£ 


-27.846 


-0.855 


2.311 


-1.843 


0.649 


-0.082 


> 3 


Rising 


100 


V€ 


-27.877 


-0.049 


0.390 


-0.138 


0.013 


0.0002 


> 3 




5000 


1/C 


-27.826 


-0.921 


2.401 


-1.902 


0.667 


-0.085 


> 3 



a Star formation history, either exponentially declining, exponentially rising, or constant. 
b Functional dependence of polynomial fit on £ = log(tage) (see Equation [T4j ■ 
c Range of ages over which the polynomial form has been fit. 
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Figure 25. Ratio of SFR to specific luminosity at 1700 A as a function of age for different star formation histories. The red curves show 
results for exponentially declining star formation histories with = 10, 50, 100, 500, and 1000 Myr. The green curve shows the result for 
a CSF history. The blue curve shows the result for an exponentially rising star formation history with t t = 100 Myr. Other rising histories 
with Tr > 100 Myr will be bracketed by the CSF and r r = 100 Myr curves. 

between L^y and S FR implies that galaxies with fainter UV luminosities also have on average lower bolometric SFRs 
(|Reddv et al.ll2OT0h . 

Second, to ensure a realistic weighting of points, we simulated a population of galaxies that has a ISchechter ( 19760 
distribution in bolometric SFR, with a ch aracteristic SFR* = 30 M B yr _1 and faint-end slope a — —1.6, based on the 
bolometric luminosity function at z ~ 2 (|Reddv fc Ste idcl 2009|). An Luv was assigned to each galaxy based on its 
bolometric SFR and the observed relationship between SFR and L\jy, as discussed above. Similarly, a stellar mass 
was assigned to each galaxy based on its SFR and the assumed intrinsic trend between SFR and M*. In assigning 
Luv an d M„, we randomly perturbed these values according to the dispersions measured in the SFR-Lyy (s» 0.5 dex) 
and SFR-M* (« 0.4 dex) relations. The final step in the simulation was to remove all galaxies with Luv ^ 1O 1O L0, 
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Figure 26. Effects of Malmquist bias on the measurements of the SFR vs. stellar mass trend, and specific SFR <f>. The solid line indicates 
an assumed intrinsic specific SFR of <j> ~ 2.4 X 10 — 9 yr — 1 . Accounting for the scatter between SFR and stellar mass results in recovered 
SFR and <f> denoted by the solid blue symbols. The hollow red symbols indicate the measured SFR and <f> when we fold in the scatter 
between SFR and UV luminosity, and impose a (spectroscopic) limit of Ltjv 
each bin of stellar mass. 



10 Lq . Errorbars reflect the dispersion in SFR and <j> in 



corresponding approximately to the spectroscopic limit at z ~ 2. The remaining galaxies are binned according their 
stellar masses. The average (or median) SFR for each bin is then calculated. 

The results of this simple simulation are shown in Figure [26l which suggests that the UV flux limit will result in 
an artificial trend of increasing specific SFR (</>) at low stellar masses, on an otherwise intrinsically flat relationship 
between </> and M*. For example, in the lowest mass bin with 9.0 < log[Af»/M Q ] < 9.5, the recovered specific SFR 
is sa 3 — 4x larger than the intrinsic value. Based on this simple simulation, and evidence that galaxies fainter than 
our spectroscopic limit have specific SFRs at M* ~ 10 9 M that are lower than what we would have inferred had we 
relied on the SED fitting of the spectroscopic sample only (see the next section), we conclude that the true median 
specific star formation is lower than that measured for the lowest mass galaxies in our sample if these SFRs are also 
determined from SED fitting (Section [7|). EEI It is evident that the magnitude of this bias increases with the intrinsic 
scatter in the SFR-M* relationship. More generally, these results highlight the importance of factoring in the UV 
luminosity limit of high-redshift dropout samples when determining mean SFRs (or specific SFRs) in bins of stellar 
mass. We show in Section [7] and Appendix [Cl how this bias can also affect our determination of the near-IR luminosity 
dependence of the mass-to-light MjLx.\p,m ratio. As such, the simulation discussed above does not tell the full story 
because we cannot directly measure the stellar masses of galaxies; the masses must be inferred from luminosity at a 
given wavelength (UV or near-IR), and therefore the position of the recovered galaxies along the stellar mass axis in 
the left panel of Figure [26] (and along both the <f> and M* axes in the right panel) may not be accurate. In Appendix IE1 
we discuss the combined effects of these biases on determinations of the SFR-M* relation. 

C. UV DEPENDENCE OF THE RELATIONSHIP BETWEEN STELLAR MASS AND NEAR-IR MAGNITUDE 

C.l. Spectroscopic Sample 

The Malmquist bias discussed in Appendix [B] can result in a slope between SFR and stellar mass that is artificially 
shallow compared to the intrinsic value. Similarly, there is a simple explanation for why this bias manifests itself 
onto the best-fit relation between stellar mass and near-IR luminosity as determined from a UV selected sample. 
Specifically, a flux-limited sample, like ours, will tend to miss galaxies with lower SFRs at a fixed (low) stellar mass 
(Section 16.41 and Appendix [B| . Therefore, the median near-IR magnitude in a given bin of stellar mass will be biased 
brighter because at a fixed stellar mass, galaxies with larger SFRs (and hence those that satisfy the UV flux limit) 
will have larger near-IR luminosities (i.e., current star formation will contribute at near-IR wavelengths). 

This effect is demonstrated clearly in Figure [27J In particular, the large scatter in near-IR luminosity (or stellar 
mass) at a given UV luminosity (see Section f7.2. 1|) . implies that there will be galaxies just below our flux limit that 
have similar stellar masses as galaxies that lie just above our flux limit. The result is that the median or average 
near-IR luminosity corresponding to a bin of stellar mass will be lower than what we would have predicted from a 
UV (bright) selected sample. Conversely, UV faint galaxies have a smaller contribution of current star formation to 
the near-IR light, and hence a given near-IR magnitude will correspond to a larger stellar mass for such galaxies. The 



26 There is also a small systematic effect in the range of specific 
SFR probed at large stellar masses. This stems from the fact that 
the mo st m assive galaxies are also more heavily star- forming (e.g., 
Figure OH) and dustier. The increase in dust obscuration means 
that at a given (large) stellar mass, we are more likely to miss 



galaxies that are heavily attenuated (because they fall below our 
UV spectroscopic limit). The effect of this bias is mitigated by 
the fact that these ultraluminous galaxies make up only a small 
fraction of the total in our sample, and thus do not significantly 
affect the fits between stellar mass and absolute magnitude. 
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Figure 27. Same as upper left panel of Figure H"3l color coded by the absolute UV magnitude (at 1700 AA), for galaxies directly detected 
in the channel 1 IRAC data. There is a clear trend between UV luminosity and the relation between stellar mass and near-IR magnitude, 
such that at a given stellar mass, UV-luminous galaxies exhibit brighter near-IR magnitudes. 

inclusion of UV fainter galaxies results in a slope between logM* and M\,\ that is less negative than the one derived 
for a sample of UV-brighter galaxies (Figure l28l This has important implications for (1) determining the l uminosity 
depend ence of the mass-to-light ratio based on a UV/optical selected samples (e.g., IGonzalez et al.l I201H ILee et all 
l2011bl ). and (2) translating the mean stellar mass found for galaxies of a given UV luminosity. 

To illustrate these points, the best-fit linear relation between log stellar mass and near-IR magnitude (log M* — M1.1) 
for galaxies in the different UV luminosity bins denoted in Figure [27] are shown in Panel A of Figure [111 Panel B 
indicates the best-fit relations when we fit for all galaxies above the faint magnitude limits of the bins in UV luminosity. 
As we add UV faint galaxies, the best-fit slope between stellar mass and near-IR magnitude becomes shallower (less 
steep). The slope of the logM* — M\,\ relation can also be affected by dust, because galaxies with larger stellar mass 
have on average larger SFRs and hence larger dust obscuration. Hence, the Panel C of Figure l28l shows that correcting 
the near-IR magnitudes for the effects of dust attenuation results in a slope of the logM* — M\,\ relation that is very 
close to 0.4, implying a dust-corrected M/L ratio at 1.1 pern that is only very weakly dependent on near-IR luminosity. 

Finally, Figure [27] sug gests that our UV selected galaxies sample approximately the full range of stellar masses 
down to Mi x — —22. [FT Therefore . Panel D of Figure l28l indicates the best-fit correlation between near-IR magnitude 
(uncorrected for dust) and stellar mass for all galaxies in our sample with Mi.i < —22 The best-fit linear relations and 
the rms dispersion of the data about these linear fits are tabulated in Table [7] The important point demonstrated in 
Figures [57] and [55] is that the systematic biases introduced by virtue of having a flux limited sample, combined with 
the effects of dust attenuation, can easily influence our interpretation of the presence of a luminosity dependence of 



2 7 No te that we have not inclu ded IR AC non-detections in Fig- 
ure 1271 As we show in Section 17.2.11 the majority (> 70%) of 
galaxies brighter than Mtjv = — 20 at z ~ 2 are detected in the 
IRAC bands. For galaxies with fainter UV luminosities, the de- 



tection fraction is lower (e.g., Figure H5). but the vast majority of 
the IRAC upper limits for such UV faint galaxies lie fainter than 
Mi.i ~ -22. 
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M AB (l.l M m) M AB (l.l M m) M~<l.l M ni) M AB (l.l M m) 

Figure 28. Panel A: Best-fit linear relations between stellar mass and near-IR magnitude for galaxies in the UV magnitude bins indicated 
in Figure [27] Panel B: Best-fit linear relations between ste llar mass and near-IR magnitude for galaxies brighter than the faint magnitude 
limit of each of the U V magnitude bins indicated in Figure 1271 Panel C: Same as Panel B, where we have corrected the near-IR magnitude 
for the affects of dust extinction, assuming the best-fit E(B — V) and either the Calzctti ct al. (2000) or SMC attenuation curve depending 
on whether the galaxy is "young" (Section [5}. Panel D: Best-fit linear relation between stellar mass and near-IR magnitude (uncorrected 
for dust) for all galaxies with Mi.i < —22.0. In all panels, the dotted lines indicate the maximum and minimum mass-to-light ratio found 
in the sample. 



Table 7 

Best-Fit Linear Relations Between Stellar Mass and Near-IR (1.1 /im) Magnitude 



Panel a 


Muv or M1.1 Range 


log[M„/M Q 


] = 


RMS 


A 


-20.5 < M UV < -19.1 
-21.0 < M uv < -20.5 
-21.5 < M uv < -21.0 
-22.0 < M uv < -21.5 
Muv < 22.0 


-(0.56±0.02)Mi 1 - 
-(0.57±0.02)Mi 1 - 
-(0.55±0.02)Mi 1 - 
-(0.56 ± 0.03)Mi.i - 
-(0.57±0.07)Mi.i - 


(2.33 ± 0.39) 
(2.58 ± 0.39) 
(2.32 ± 0.46) 
(2.51 ± 0.67) 
(2.80 ± 1.64) 


0.19 
0.17 
0.16 
0.16 
0.19 


B 


A%v < -19.1 
M uv < -20.5 
M uv < -21.0 
M uv < -21.5 


-(0.52 ± 0.01)Mi 1 - 
-(0.53±0.01)Mi 1 - 
-(0.54±0.02)Mi 1 - 
-(0.55 ± 0.03)Mi.i - 


(1.49 ±0.24) 
(1.75 ± 0.26) 
(1.98 ±0.35) 
(2.36 ±0.62) 


0.20 
0.18 
0.17 
0.17 


C 


M uv < -19.1 
M uv < -20.5 
M uv < -21.0 
M uv < -21.5 
M uv < -22.0 


-(0.46±0.01)Mf° r - 
-(0.46 ± O.OlWfJ - 
-(0.46±0.02)Mf ) 1 r - 

-(0.45 ± om)Mi°i - 

-(0.43 ± 0.07^^ - 


(0.15 ±0.23) 
(0.26 ±0.27) 
(0.37 ±0.36) 
(0.32 ± 0.64) 
(0.31 ± 1.67) 


0.22 
0.20 
0.19 
0.20 
0.24 


D 


M1.1 < -22.0 


-(0.56 ± 0.09)Mi.i - 


(2.42 ± 1.94) 


0.12 



a Indicates panel of Figure [28] that shows the best- fit linear relations between M, 
and M1.1. 



the M/L ratio. 

C.2. UV Faint Sample 

We have several options for inferring the stellar masses of UV faint galaxies that lie below our spectroscopic limit. 
One possibility is to simple assume the same relationship between near-IR magnitude and stellar mass as was found for 
the faintest UV luminosity bin in our spectroscopic sample. A second option is to attempt to fit for the average SED 
of these UV faint galaxies. This option is beyond the scope of this paper and is not explored further. A third option is 
to extrapolate UV luminosity dependence of the log M* — M1.1 relation to UV faint galaxies. While the marginalized 
uncertainties in the slope and intercept of the logM* — Mi. 1,0.9 relations for the different UV luminosity bins in the 
spectroscopic sample suggest very little evolution (Table [7]) , there is some indication that the relation systematically 
shifts in normalization towards higher masses at fainter UV- luminosities (Panel A of Figure [28|) . as would be expected 
if these UV faint galaxies have a lower contribution of current star formation to the near-IR magnitude (sec discussion 
above and in Section l7.2.1l) . By extrapolation from the UV-bright bins, the intercepts of the logM* — Mi. 1,0.9 relation 
are computed as b = 0.14Mtjv + 0.55. The slope is kept fixed at S log[M*/M ]/<5Mi. 1,0.9 = —0.56 (similar to that 
found for the more UV luminous bins). Converting the stacked IRAC channel 1 magnitudes listed in Tabled we 
find that the median stellar masses obtained in this way are log[M*/M Q ] = 9.02 and 8.71 for the two UV faint bins 
at z ~ 2. Similarly, we find ^[M^/M©] = 9.19 and 9.12 for the two faint bins at z ~ 3. The uncertainties in 
these stellar masses include: (1) measurement uncertainties in the stacked IRAC fluxes, (2) statistical uncertainties in 
the fit between near-IR magnitude and stellar mass, and (3) systematic uncertainties in the relation used to convert 
between near-IR magnitude and stellar mass (e.g., arising from the use of other of the other options mentioned above). 
Combining the errors from these effects in quadrature results in stellar mass uncertainties of rs 0.15 — 0.29 dex. 
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Table 8 

Stacked IRAC Magnitudes for Spectroscopic Sample 



Redshift Interval M(1700A) Range 3.6 /im a 4.5 ^tm a 5.8 /im" 8.0Aim a 



1.4 < 


z < 2.7 


-22 


.5 


-22 


.0 


-23 


.92 ± 0. 


13 


(4) 










-24.06 ±0.16 


(2) 






-22 


.0 


-21 


.5 


-22 


.87 ±0. 


21 


(28) 


-22.93 ±0.19 


(26) 


-23.19 ±0.20 


(12) 


-22.95 ±0.20 


(22) 






-21 


.5 


-21 


.0 


-22 


.20 ±0. 


17 


(72) 


-22.27 ±0.17 


(48) 


-22.41 ±0.19 


(31) 


-22.25 ±0.21 


(55) 






-21 


.0 


-20 


.5 


-21 


.95 ± 0. 


20 


(192) 


-21.85 ±0.19 


(138) 


-22.14 ±0.20 


(78) 


-21.78 ±0.23 


(147) 






-20 


.5 


-20 


.0 


-21 


.61 ± 0. 


19 


(199) 


-21.79 ±0.20 


(180) 


-21.99 ±0.21 


(93) 


-21.68 ±0.22 


(189) 






-20 


.0 


-19 


.5 


-21 


.14 ±0. 


21 


(113) 


-21.15 ±0.19 


(79) 


-21.35 ±0.26 


(39) 


-21.02 ±0.29 


(93) 






-19 


.5 


-19 


.0 


-20 


.85 ± 0. 


21 


(22) 


-20.96 ±0.26 


(22) 










2.7 < 


z < 3.7 


-23 


.0 


-22 


.5 


-23 


.51 ± 0. 


27 


(2) 


















-22 


.5 


-22 


.0 


-22 


.45 ± 0. 


13 


(13) 


-23.18 ±0.12 


(6) 


-23.06 ±0.32 


(6) 


-23.54 ±0.21 


(6) 






-22 


.0 


-21 


.5 


-22 


.35 ± 0. 


12 


(45) 


-22.49 ±0.13 


(34) 


-22.32 ±0.25 


(21) 


-22.57 ±0.16 


(34) 






-21 


.5 


-21 


.0 


-21 


.83 ± 0. 


13 


(104) 


-21.83 ±0.12 


(78) 


-21.87 ±0.18 


(52) 


-21.99 ±0.19 


(79) 






-21 


.0 


-20 


.5 


-21 


.37 ±0. 


12 


(108) 


-21.49 ±0.13 


(84) 


-21.48 ±0.25 


(65) 


-21.13 ±0.30 


(81) 






-20 


.5 


-20 


.0 


-20 


.97 ±0. 


16 


(68) 


-21.12 ±0.12 


(57) 











a Uncertainties in absolute magnitude reflect the stacked flux measurement uncertainty combined in quadrature with the dispersion 
in absolute magnitude given the range of redshifts of objects in each bin. Parentheses indicate the number of galaxies in the 
stack. 



D. IRAC STACKING METHOD 

As noted in Section [3J we modeled the stellar populations only for those galaxies that were directly detected at 
wavelengths longward of rest- frame 4000 A to ensure that robust mass determinations included information from the 
age-sensitive Balmer and 4000 A breaks. However, a substantial fraction of galaxies in the spectroscopic sample are 
undetected at these longer wavelengths, particularly those with fainter UV luminosities, as demonstrated in Figure [T51 
For an unbiased view of the mass distribution, we stacked the IRAC data in bins of UV luminosity. 

Our stacking analysis proceeded as follows. Cutouts around each galaxy in our sample were assembled from the 
IRAC images. We included in the stack any galaxy that (1) does not show evidence of AGN, (2) has a spectroscopic 
redshift in the range 1.4 < z < 3.7, (3) has an IRAC exposure time of more than 7200 sec, (4) is not blended with any 
nearby neighbor, and (5) lies more than 2'.'4 away from any nearby source as identified in higher resolution optical and 
near-IR data. The cutouts were then median combined, though average combination yielded similar results. Given 
the field-dependent variation of the IRAC PSF, we constructed a "stacked PSF" made in exactly the same way as 
the regular stack — i.e., for each galaxy that went into the stack, we also stacked the PSF appropriate for the field 
in which that galaxy lies. Photometry on the stacked image was measured by fitting the stacked PSF to the stacked 
IRAC signal. Table [8] summarizes the stacked photometry, converted to absolute magnitudes based on the median 
redshifts of the objects that went into the stacks. 

E. COMBINED EFFECT OF BIASES ON THE SFR-Af» RELATION 

In this section, we discuss the combined affect of the biases discussed in Appendices [B] and [C] To recap, there is 
Malmquist bias of selecting galaxies with larger star formation rates at a given stellar mass in a flux limited sample. 
Second, there is a bias introduced by the conversion that one used to convert near-IR luminosity to stellar mass. 
To investigate jointly these effects, we built upon the simulations discussed in Appendix [Bj For each galaxy in the 
simulation, we assigned it a near-IR magnitude M\.\ based on its UV magnitude, using the correlation between these 
two quantities, taking into account upper limits in near-IR magnitude with the EM parametric estimator (Figure [T4"l) : 

Mi.i = 0.90Afi7oo - 3.41 (16) 

with a dispersion of ldex. This near-IR magnitude was then converted to a "measured" stellar mass, M™ cas , using 
three different options. In the first option, M™ oas is estimated for each galaxy using the logM* — Mn relation 
appropriate for the UV luminosity of that galaxy (Panel A of Figure EH and Table [7]) , and assuming that galaxies 
fainter than M1700 = —19.1 follow the relation for galaxies with —20.5 < M1700 < —19.1 . In the second option, M™ cas 
is estimated for each galaxy using the log M* — M\.\ relation found for all galaxies in our spectroscopic sample with 
Mnoo < —19.1 (Panel B of Figure [28l and Table [7]). In the third option, M™ eas is estimated for each galaxy using a 
logAf* — M1.1 relation with a fixed slope of (51og[Af*/M Q ]/(5Mi.i = —0.56 and an intercept that evolves with M1700 
(see Appendix ICj). 

The comparison between the simulated (M^ lm ) and measured stellar masses is shown in Panel (a) of Figure [29j 
For the first two options, there is a tendency to underestimate Mf m because we attribute too large a fraction of the 
near-IR light to star formation rather than stellar mass. With the third option, we obtain a reasonable agreement 
between M^ lm and M™ cas . Applying a magnitude cut results in an artificial offset between Af™ eas and M* lm because 
fainter galaxies will be excluded by the magnitude limit, leaving only brighter (and more massive simulated) galaxies 
in the bins of M™ cas . The middle row shows the recovered SFR-M* relation, where galaxies are binned by M™ cas . 
Because options 1 and 2 result in an underestimation of stellar mass, there is tendency to predict a shallower slope of 
the SFR-M* relation. Introducing a 1Z < 27.0 limit then causes us to overestimate the median SFR in bins of M™ cas , 
due to Malmquist bias (Appendix [Bj . Finally, the bottom row shows the results when binning by SFR, instead of 



34 



10 



10 



10 9 



E 

«* 

58 10 8 





03 
CO 



03 

Cl- 
og 



10 7 ^ 



10 1 r 



10° r 



10" 1 - 



S 10° 



10 



-1 =- 



10 



-2 



i i 1 1 niij — i i 1 1 niij — r~r 

(a) No Mag 




10-2 r+444 



10 1 r 



10 7 



10 8 10 9 



10 7 



^rneas 



(M ) 



10 8 
Mmeas 



i 

10 9 10 10 



(M ) 



Figure 29. Comparison between simulated (black lines) and measured stellar mass (Af™ cas ), binned by the latter, for galaxies in the full 
simulated sample (top row, panel a) and those brighter than TZ = 27.0 (top row, panel b). The cyan, orange, and red lines assume options 
one, two, and three described in the text for converting near-IR, magnitude to stellar mass. The middle row indicates the same for the 
median input SFR in bins of AfJ ncas . The bottom row indicates the same, except binned in terms of SFR. 

Mf eas . In this case, as expected, the Malmquist bias is less noticeable, because at a given low SFR (or faint UV 
luminosity) we can probe the full range of stellar mass, or at least quantify the average stellar mass using stacking 
analyses (Figure [U and Appendix [D]) . The results from the simulations underscore how the obvious Malmquist bias 
and the subtler bias arising from the conversion between near-IR light and stellar mass can affect our interpretation 
of the SFR-Af* relation. In practice, estimating the SFR-Af, relation in bins of SFR may yield more accurate results 
that are largely immune to the effects of Malmquist bias, though one will still have to account for systematics in the 
conversion between near-IR light and stellar mass, particularly if SED-fitting for all the objects in a sample is not an 
option. In our case, for the UV faint samples, we determined median stellar masses in bins of UV luminosity, which 
roughly translates into bins of SFR, so our estimates should be largely unaffected by the Malmquist bias discussed 
above. 



